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A BSTR ACT 

Predicting the thermal structure of the oceans is of 
importance to the Naval tactician, logistician, or search 
and rescue coordinator. Understanding the structure of the 
oceans provides valuable insights to those who must utilize 
the oceanic environment effectively in their day to day 
operations . 

Today, recent information about an area is limited to 
point observations of single bathythermographs. Few models 
produce an accurate picture of the ocean environment that 
can be used for updating tactics to conform to a changing 
situation. Producing a reliable prediction of conditions 
for a large area, while using limited resources, is the 
basic objective of this paper. 

Satellite infrared imaging of the ocean surface has been 
used effectively to map sea surface t emperat ure patterns. 
Such sea surface temperature patterns can be used, along 
with climatology, to identify subsurface thermal structure 
in an ocean area according to results of this study. More 
accurate inputs can be made to range dependent acoustic 
prediction models, thus improving ‘he ASW environmental 
predictions available to fleet users. 
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I. INTRODUCTION 

A. PREDICTION OF OCEAN THERMAL STRUCTURE 

Attempts to predict the vertical temperature profile in 
a large body of water can be quite varied. A climatological 
atlas may be used and a profile chosen which corresponds to 
the location and time of interest. The climatology may be 
updated using an interpolation from one climatological 
period to another. A recent measurement in an area can be 
used to modify the climatological prediction, thus producing 
a more time-sensitive prediction. Measurements which may be 
used to to modify climatological predictions are temperature 
structure in the upper ocean and sea surface temperature 
(SST) . 

The bathythermograph (3T) provides a temperature profile 
at a point, while satellite measurements can provide recent 
sea surface temperature patterns over an area of interest. 
A combination of climatology, recently-measured temperature 
structure (ET) , and satellite-derived SST could provide a 
very reasonable prediction of the vertical temperature 
profile for a body of water . 

3. PROBLEM OF ACOUSTIC PREDICTIONS 

Accurate predictions cf transmission loss depend on 
accurate specifications of wafer masses along an acoustic 
path. The primary physically-measurable quantity which 
affects sound speed (and thus transmission loss) is tempera- 
ture. Temperature structure of the water is an important 
factor in acoustic losses along a transmission path. There 
is no way to gather temperature profiles rapidly along a 
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desired track. To measure temperatures along the tracks 
analyzed for the experiment discussed in this thesis 
required a nine-hour mission by a dedicated patrol aircraft. 
The operational fleet user does not have the luxury of dedi- 
cated air assets to provide a detailed thermal description 
of an area of acoustic interest. 

Obviously, the employment of a dedicated patrol aircraft 
loaded with a store of expendable bathythermographs is not a 
realistic solution to the problem of predicting acoustic 
properties in various locations of the world. The use of 
numerical models, relying on climatological files, and modi- 
fied by recent observations, provides a first guess in 
predicting thermal structure. However, these models require 
time to produce results, and may therefore not accurately 
present the actual thermal structure. 

If a coupling exists between conditions at the air-sea 
boundary and the depth of the mixed layer, satellites can 
provide a means of rapidly surveying the sea surface in an 
oceanic operating area. The satellite sea surface 
observations could then be used to estimate the vertical 
temperature structure throughout the operating area. If a 
good linkage exists between sea surface conditions and the 
mixed layer depth, the accuracy and timeliness of acoustic 
predictions could be improved. Such ocean surface 
conditions observable remotely by satellite include sea 
surface temperature, ocean color, and topography. 

C. EXPERIMENTAL BASIS FOR THIS THESIS 

In November and December of 1980, the Naval Postgraduate 
School conducted the Acoustic Storm Transfer and Response 
Experiment (ASTRSX) . This experiment was conducted as part 
of the first Storm Transfer and Response Experiment (STRSX) , 
a joint investigation by United States and Canada to 
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determine the effects on the air-sea boundary of North 
Pacific storms, and the converse [Ref. 1]. Figure 1.1 shows 
the location of the A SIR EX area in the context of the 
eastern North Pacific Ocean. Three flight tracks were 

examined within this area. Each track covered a diagonal 
through the central area of the square from the southeast to 
the northwest. 
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II. CLIMATOLOGY OF THE EASTERN NORTH PACIFIC OCEAN 

A. OVERVIEW OF THE EASTERN NORTH PACIFIC OCEAN 

Much of the research conducted in the Eastern North 
Pacific Ocean relies on data collected at Ocean Weather 
Station "P" (OWS-P) , located at 50N, 145W, in the center of 
the Gulf of Alaska. The research project reported here 
operated in an area of the North Pacific to the east and 
south of OWS-P. Tabata [Ref. 2] concluded that the charac- 
teristics of OWS-P were representative of the wafers in the 
northeast Pacific region 

Ocean Station "P” is situated north of the Subarctic 
Current, well to the northwest of the divergence which is 
formed when the West Wind Drift approaches the North 
American continent [Ref. 3]. The California Current is the 
southerly component of the West Wind Drift after it divides, 
while the northerly component forms the Alaska Current and 
continues i-s flew into the Alaska Gyre [Ref. 4]. 

The climat ologic al description of "-he project area 
depends on a combination of multiple-year observations from 
OWS-P and observations within the area itself. The region 
from which -he data were taken includes the subarctic front 
and the source of the California Current. The California 
Current is a permanent feature in the observation area. 

The observation area may be described by three regions: 
the Coastal, or Transition Region, extends from the coast 
out to 130W; the Subarctic Watermass extends seaward beyond 
the coastal zone to the north of 42N; and the Central 
Pacific Watermass is south of 4 2N. [Ref. 5 1. Boundaries 
separating these zones are dynamic, and demarcation lines 
serve only as initial positions for observing variability. 
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Figure 2.1 Watermasses of the Eastern North Pacific Ocean 
[Ref. 6 j . 
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A seasonal change in the location of these boundaries 
can be noted from observations of climatological data. 
These seasonal changes can be related, in part, to weather 
conditions in the Northeast Pacific. In the winner, the 
Aleutian Low dominates the atmospheric sea level pressure 
patterns, and in the summer the region is under the 
influence of the North Pacific High. These systems both 
drive the prevailing westerly winds, which give rise to the 
West Wind Drift in the upper ocean [Ref. 7]. 

Although the boundaries do vary, the variation of the 
boundary between the Subarctic and the Subtropic watermasses 
is low. The position of the front was "relatively 
constant", between 40N and 42N, as observed by Dodimsad et 
al . [Ref. 8]. Roden [Ref. 9] found the position of the 
front to vary from 40 N to 45N, and stated that this frontal 
position was highly dependent on the prevailing wind fields. 
This seems reasonable for a latitudinal boundary as the wind 
direction remains predominately westerly from month to 
month, while wind speed in winter is generally twice that of 
summer [Ref. 10 ]. 

Dodimead £t ai- found than the coastal zone was much 
smaller in winter than in summer [Ref. 11]- This could be 
due in part to an increase in wind speed in winner (compared 
to summer) along a longi tud indai boundary that is termir.a x ed 
by the coastal zone and North America. 

B. WATER MASS CHARACTERISTICS 

The region used for A5TREX contains three water masses. 
These are: 1) Coastal Water extending from the coast to the 
California Current (130W); 2) Eastern Subarctic Pacific 
Water extending southward to about 42N; and 3) Subtropic 
Pacific Water to the south of 42N [Ref- 12]- 
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Dodimead et al. reviewed the area in terms of various 
domains. This description shows three zones in + he vertical 
salinity structure: the upper zone, located above a perma- 
nent halocline, displays a seasonal variation; the permanent 
halocline (100-200m) represents the second zone whose 
vertical extent varies with latitude; below the halocline is 
the lower zone, which begins at the salinity surface (S = 



33. 8 g/kg) , and has 


gradually-increasing salinity 


and 


decreasing temperature 


with depth. 


It displays 


less 


seasonal variation in 


temperature and 


salinity than 


the 


upper zone. [Ref. 13]. 








In the upper zone. 


the region is 


influenced by 


both 


Coastal and Transitiona 


1 Domains. The 


lower zone ca 


n be 


partitioned into the Central Subarctic, 


Transitional 


and 



California Undercurrent Domains [Ref. 14]. 

The Tabata and Dodimead et al. descriptions vary in 
approach, but results are similar for the vatermass 
structure of the ASTREX region. 

1 * Wat er mass A na ly sis 

The Coastal Watermass, referred tc as the Transition 
Region by Sverdrup, Johnson and Fleming (1942), extends from 
the coast of California to about 130w Longitude, where it is 
bounded by the Eastern North Pacific Central Water (or 
Subtropic Pacific Water according to Tully, 1964). 

Both the southward California Current and the 
inshore countercurrent (Davidson Current) make the 
Transition Zone highly variable. By contras: , booh the 
Subarctic Pacific Water and the Subtropic Pacific Water have 
much less variability over time. 

Predominately during spring and summer, upwelling 
occurs off the California coast. In regions of upwelling, 
the sea surface temperatures in summer may be lower than in 
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winter. The Davidson countercurrent continues to flow 
northward during the upwelling season, but at depths greater 
than 200 meters. During the winter, the Davidson Current 
may extend to the surface and flow northward no about 48N. 

An examination of monthly mean depths to the ^cp of 
the thermocline shows a depth variation of 15 to 60 meters 
at the end cf the upwelling season ( Jul y-September ) . A vari- 
ation of 45 to 90 meters in the depth of the thermocline can 
be noted at the end of the winter season (January-February) . 
Additionally, in all seasons, the mixed layer deepens with 
increasing distance from the coast. [Ref. 16], 

Close to the coast, in the Transition Zone, salinity 
and temperature structures will be affected by both 
upwelling and ccastal run-off. With increasing distance 
from the coast, a more stable structure develops. West of 
130W, salinities at the surface are 32.5 to 33.0 g/kg, 
increasing to about 33.0 to 33.5 g/kg at 120 meters depth. 

The Eastern Subarctic Pacific Watermass is charac- 
terized by three layers. The upper layer (surface to 
approximately 100 meters) is seasonal in structure. In 
winter this layer is well-mixed, but in summer a sharp ther- 



mccline and a halocline form in this layer. 



The second 



layer is the "princ ip al" halocline (from 100 to 200 meters) 
which marks the transition from the seasonal layer to the 
"stable" deep layer. In the deep or lower layer (depths 
greater than 200 meters) , the salinity increases with depth 
and the temperature decreases. At 1000 meters depth, the 
temperature is 2.8C and salinity is 34.4 g/kg. 

By contrast, the Subtropic Pacific Water has no 
permanent halocline. Rather, it shows a salinity decreasing 
with depth to a minimum between 200 and 300 meters. 
Temperature is the dominating factor in the formation of the 
seasonal layer. The upper zone, which displays seasonal 
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temperature variations, extends to a depth of aoout 150 
meters, beyond which a permanent thermccline exists. This 
seasonal zone is isothermal to a dep^h of 150 meters in 
early spring, shallowing to a depth of 40 to 60 meters 
during the summer and into fall. 

2 • Do main Ana ly s is 

By partitioning the area into well-defined domains, 
a more detailed description of the wafer masses can be made. 
However, there are few clear-cut boundaries in either the 
upper or lower zones, and the halocline vanishes in certain 
situations [Ref- 17]. 

The upper zone of the test region is a combination 
of Transitional and Coastal Domains. The Coastal Domain has 
a thermal structure which is very dependent on location. 
Such location-dependent processes as river run-off, heating 
and cooling, and wind cause the Coastal Domain to be highly 
variable from place to place [Ref. 19]. 

The Transitional Domain can be identified by rela- 
tively warm waters (T>15c in summer and T>7c in winter) , and 
relatively high salinities (S>32.2 g/kg at the surface and 
S>33.4 g/kg at the top of the halocline). In the test 
region however, the salinities may be somewhat lower 
[Ref. 19] (as discussed earlier, mean salinities at the 
surface in the test region were 32.5 to 33.0 g/kg; at the 
top of the halocline the salinities were 33.0 to 33.5 g/kg) . 

Based on the domain description, the upper zone of 
the test region can be characteriz-d by a layer cf water 
approximately one hundred meters deep, which varies 
seasonally. This upper zone shows an increase of salinity 
with depth to the halocline when it exists. An even 
stronger increase cf salinity occurs through the halocline 
(.003 g/kg/m < G3AD ( S) < .008 g/kg/m). At the bottom cf 
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125°W 




Figure 2.3 Domains in the ASTRSX Region. 
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Figure 2.4 



Upper Zone Domains 



[Ref. 20]. 
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the halo cline the salinity approaches 33.8 g/kg. Also, a 
temperature minimum may occur at the bottom of the upper 
layer [Ref- 21]. 

The upper zone responds to atmospheric affects and 
shows much more variability than does the lower zone. The 
halocline acts as a ’’permanent boundary" between the upper 
zone and the lower zone. Compared to the upper zone, the 
lower zone has a stable salinity structure. [Ref- 22]- 

A second halccline forms in the upper zone during 
the summer. This halocline begins at 10 J -o 30 meters depth 
and is due to the mixing effects of the light summer winds. 
Below this secondary halocline (which can extend to 75 
meters), the water is nearly isohaline down to the principal 
halocline (100 to 20 0 meters dep-^h in the Subarctic Water) 



[Ref- 23]. 

The lower zone of the test region can be partitioned 
into (at most) three domains, depending on the season. 
These lower domains are the Central Subarctic, Transitional, 
and the California Undercurrent [Ref. 24]- 

The Central Subarctic Domain is typical of the lower 
layer water in the northwest corner of the test region. 
This is also the water type associated with OWS-P. Dodimead 



et a 1 . identify this domain as water underlying a 
well-defined halocline. At OWS-P, the halocline is defined 
by a salinity gradient of 0.91 g/kg per meter. The salinity 
increases with depth very slowly in the lower layer 
[Ref. 25]. 



The Transitional Domain underlies 



its 



sur race 



counterpart throughout the southern half of the test region 
and also extends through the northeast quarter of the 
region. Near the coast, both this domain and the California 
Undercurrent can exist. The Transitional Domain is strongly 
influenced by the major geostrophic currents in the area. 
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Figure 2.5 



.rovth/decl me of the Thermocline at OWS-P 
Ref. 26] . 
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Figure 2.6 Lower Zone Domains [Ref. 27] 
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It is a narrow band underlying the West Wind Drift and 
broadens as this current separates tc form the Alaskan and 
California currents. In the Transitional Domain the 
halocline is less pronounced than in the Central Subarctic 
Domain. The Transitional Domain bridges the gap between the 
strongly-layered Subarctic Water and the Subtropical Water 
which has no distinct halocline [Ref. 28]. 

The third lower domain present in the test ar^a is 
the California Undercurrent Domain. As the name implies, 
this domain is influenced principally by the subsurface 
countercurrent which runs northward near the coast. The 
countercurrent originates in the warm water of the Eastern 
Central Pacific Water. Dodimead et al . characterizes 
California Undercurrent water by temperatures greater than 
6C on the 34 g/kg surface of salinity, as opposed to the 
Transitional Domain defined at a salinity of 33.8 g/kg with 
temperatures less than 6.5C [Ref. 29]. 

3 . Te m pera ture S tru ct u r a 

The salinity structure of the eastern North Pacific 
provides a relatively stable basis for defining water masses 
in the area. However, characterization of water masses 
requires a salinity-temperature relationship as well. The 
true variability of acoustic conditions in the eastern North 
Pacific must come from periodic and/or random variations in 
temperature which are superimposed on the basic density 
structure determined by salinity. Of course, salinity 
variations may be of importance in the coastal regions too, 
where siainificant intrusions of fresh wafer or more saline 
upwelled wafer may occur. 

At OWS-P, the temperature of the upper zone has an 
annual periodic variation, reaching a maximum in the late 
summer, and a minimum in late winter. The maximum annual 
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temperature variation noted between 1956 and 1962 was about 
9C. Variations of temperature in the entire upper zone 



followed the same 


cycle 


as 


sea surface 


temperature. 


appearing i o be in 


phase 


with 


the periodic 


sea surface 


temperature variation 


[Ref. 


30]. 






As with sali 


n it v , 


the 


temperature var 


iations are 



greatest in the upper zone. The temperature range decreases 
from a maximum of 9C in the upper 59 meters to less than 3C 
variation at 125 meters depth. Below the haloclir.e, the 
small temperature variation with depth (gradient) should be 
reasonably well predicted from monthly climatological data. 

In the upper zone, a seasonal thermocline begins to 
develop in spring with the advent of the heating season. 
The isothermal March layer becomes a two-layer medium 
through April when a seasonal thermocline begins to develop 
at abcu* 30 meters. This seasonal thermocline reaches a 
maximum depth of about 50 meters, depending on the strength 
of the prevailing winds, but can be expected to occur 
between 30 and 50 meters during the spring and summer. The 
thermocline shows the strongest gradient in late summer 
(Au gust-Se ptember) , when it is at its shallowest depth (30 
meters) [ Ref. 31]. 

With the advent of autumn, the cooling season 
begins. The seasonal thermocline deepens throughout the 
cooling season due to the effects of colder surface water 
and increased mixing caused by the strong winter winds. 

The seasonal thermocline forms, coincident with the 
formation of the secondary or seasoral haloclir.e. This 
thermocline is very stable, and hence, the structure of the 
upper zone is quite predictable [Ref. 32]. The seasonal 
thermocline and permanent haloclir.e locations can be used as 
inputs to predict the vertical structure of the upper 100 
meter layer during the heating and cooling seasons. Sea 
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surface temperatures may provide input regarding the 
strength of the seasonal thermocline (gradient) since both 
the layer thickness and the temperature at the bottom of the 
thermocline remain relatively constant. 

Characteristics of the Subarctic Wafer are found in 
the waters which comprise the Transitional Domain. Wafers 
in the Coastal Domain, on the other hand, will vary from the 
Subarctic Water structure. Coastal Waters remain very near 
the coast in winter, but extend further seaward in the 
summer. The primary effect of mingling Coastal Waters with 
Transition Water will be on salinity structure. The Coastal 
Water tends to be less saline than Transition Water. A 
winter section of temperature and salinity taken through the 
test region in 1959 indicated a narrow zone of coastal 
influence (to about 130W) , which showed as an increase in 
temperature and decrease in salinity from the Transitional 
Domain. A summer section at the same latitude showed a 
strong coastal influence out no 130W; the temperatures were 
much warmer and salinities much lower in the Coastal Domain 
than in the Transitional Domain. 

C. METEOROLOGY OF THE EASTERN NORTH PACIFIC OCEAN 

The discussion of the meteorology of the ASTRSX region 
and the northeast Pacific Ocean will be limited to factors 
which affect spaceborne instruments, rather than those 
factors which influence the circulation and heat transfer in 
the area. 

The most obvious path interference encountered by a 
spaceborne instrument is atmospheric moisture in the form of 
clouds. During ASTREX, clouds obscured most of the satel- 
lite field of view on every measurement day. Since the 
purpose of the experiment involved studying the effects of 
North Pacific storms, cloudiness was an unavoidable obstacle 
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for the secondary mission of comparing satellite data to In 
Situ observations. Winter cloudiness in the North Pacific 
is a major problem for visual and infrared imaging of the 
ocean surface. 

The meteorological phenomena which directly affec 1 * 
space-based ocean viewing instruments are aerosols in the 
form of moisture and salts. Moisture, in its obvious forms 
of clouds and rain, severely degrades passive observation 
instruments which rely on visual or infrared energy. 
Passive microwave radiometers suffer degradation when the 
viewing area includes severe rains. 

The Navy Marine Climatic Atlas of the World [Ref. 33] 
was consulted to determine climatic conditions in the 
Eastern North Pacific in November and December. For both 
months, cloudiness greater than 4/8 can be expected about 
four days in five. This represents a severe degradation for 
visual and infrared imaging of the North Pacific waters in 
November and December. 

Precipitation reported in the Atlas increases from 19 
percent in November to about 24 percent in December. In 
other words, rain occurs about one day in five in November 
and one day in four in December. Rain constitutes an inter- 
ference effect for passive microwave systems and total 
degradation for visual and infrared systems. 

Relative humidity was reported as greater than 90 
percent one day in five in November and one day in four in 
December. This follows the precipitation variation. Maul 
and Sidran [Ref. 34] reported on effects of water vapor at 
45N and 60N. From their report, wafer vapor content in the 
atmosphere ever the Eastern North Pacific (between 4Q-50N) 
can cause a temperature error of the order of 2 degrees 
Kelvin in winter and 6 degrees Kelvin in summer. These 
errors represent the difference between unconnected 
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satellite sea surface temperature observations and 
sea surface temperatures. 



act ual 



33 



Ill 



. SATELLITES IN OCEANOGRAPHY 

A. INSTRUMENTS FOR OBSERVING THE OCEANS 
1 • Visual and In fra rad Imaging 

Several series of satellites have both visual and 
infrared capabilities. Among these satellites are the GOES 
series, NIMBUS series, and NOAA series. Stewart (1981) 
presents an excellent summary of satellite instruments used 
for oceanography. He also lists the satellites which carry 
oceanographic instruments and the specific instruments each 
satellite carries [Ref. 35]. 

The Geostationary Observational Environmental 
Satellite (GOES) system is composed of two satellites in 
geosynchronous earth orbits over the equator. The satel- 
lites, designated GOES-East (subpoint at 75 degrees west 
Longitude) and GOES-West (subpoint at 135 degrees west 
Longitude) , scan the earth, with overlapping coverage over 
the Americas, from Africa to the western Pacific Ocean 
[Ref. 36]. 

The GOES satellites, using the visible and infrared 
spin scan radiometer (VISSR), can scan the full disk of the 
earth in 18.2 minutes and generate an image every thirty 
minutes [Ref. 37]. Figure 3.1 is a visual image produced by 
the GOES VISSR. Table I lists the spectral channels and 
resolutions for the VISSR. Figure 3. 2 presents the atmos- 
pheric transmission windows for the visible and infrared 
portions of the electromagnetic spectrum. 

GOES can provide near-real-time coverage of oceanic 
areas in the Atlantic and Pacific Oceans. GOES imagery was 
utilized during the STREX experiment as a meteorological 
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T ABLE I 



Visible and Infrared Spin Scan Radiometer (VISSR) 

Chann elization. 



Channel Spatial Resolution 
{at subpoint) 



Wavelength 
(micr om eters) 



Vis ib le 
Infrared 



1 km 
8 km 



0.55 - 0.75 
10.5 - 12.6 



Note: 

1. Fett gt 
channel 
for the 



al . list spatial resolution for the visible 
as 0.5 nautical miles, and 5 nautical miles 
infrared channel. 



aid. GOSS images of sea surface temperature are useful for 
studying oceanic circulation, but the spatial resolution is 
not as good (8 km) as the resolution from satellites in low 
ear- 1 - h orbits, such as NOAA-6 (1.1 km resolution at infrared 
wavelengths) . 

The Nimbus-7 spacecraft carries the Coastal Zone 
Color Scanner (CZCS) , which is a visual radiation instrument 
centered in specific wavelength bands to measure 
ccncentra- ions of phytoplankton in the ocean. Table II 
shows the wavelengths for indicated CZCS channels. 

The Advanced Very High Resolution Radiometer (AVHRR) 
is a five-channel instrument used on the NOAA-o and NOAA-7 
spacecraft. The AVHRR on the NOAA-6 has one visual channel 
and three infrared channels (the fifth channel is a repeat 
of channel four). NOAA-7 utilizes all five channels. The 
AVHRR has a ground resolution of 1.1 km in all channels. In 
the daytime, the first two channels (visible and far red to 
near IR) can be used to filter cloud-contaminated infrared 
measurements. Channel three (3.5-3.93 urn) can only be used 
at night because of severe contamination by reflected solar 
energy. Channel four (1 0. 5 - 1 1 . 5um) can be used day or night 
for measuring sea surface temperature. 
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Figure 3. 1 



GOES Visible Image 



(Full Disk) 



[Ref. 38], 
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*»VELENCrH - ki m 



Figure 3.2 Atmospheric Tr ansmittance (Visible and Infrared 
Spectra) [Ref. 39], 
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TABLE II 



Coastal Zone Color Scanner (CZCS) Channelization. 



Chann el 


Wavelength (mi 


cr oa e ters) 


1 


0. 433 - 0.453 


(blue) 


2 


0. 51 0 - 0.530 


(green) 


3 


0. 540 - 0.560 


(ar sen) 


4 


0. 660 - 0.680 


(red) 


5 


0. 700 - 0. 800 


(far rad/ 


6 


1 0. 5 - 12. 5 


(far IR) 



IR) 



Not e: 

1. Channel 6 failed in 1979. 

2. Resolution is approximately 1 km for all channels. 



Table III is adapted from booh Reference 36 and 37. 
It lists the AVHRR spectral channels for the AVHRR/1 and 
AVHER/2. 



TAELS III 

Advanced Very High Resolution Radiometer (AVHRR) 

Chan r.elization 



Channel Wavelength (micrometers) 



1 

2 

3 

4 

5 



0.58 - 0. 68 
0.725 - 1.10 
3.55 - 3. 93 

10.5 - 11.5 

11.5 - 12.5 



i visible) 
far red-near IR) 
middle I R) 
far IR) 
far IS) 



Not e : 

1. Spatial resolution at satellite subpoint is 
0^5 nautical miles (1.1 km) . 

2. Channel 4 wavelencrth is 10.3 - 11.3 um for all 
AVHRR/2 instruments. 

3. Channel 5 added to further enhance sea surface 
measurements in the tropics. 
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The noise equivalent temperature (NETD) of the AVHRR 
is less than 0.12 degrees Kelvin for a 300 degree Kelvin 
scene [Ref. 40]. The temperature sensitivity for the chan- 
nels is approximately 0. 5 -degr ees-Kel vin when viewing a 
300-degr ee-Kelvin scene. Since the noise equivalent temper- 
ature of the system is less than the temperature sensitivity 
of the channels, identification of thermal gradients on th® 
order of 0 . 5-degrees- Kelvin per kilometer is feasible. 

In Chapter two, the meteorological phenomena which 
interfere with visible and infrared imagery of the ocean 
surface were discussed. The following discussion deals with 
the AVHRR and techniques used to compensate for atmospheric 
effects. The AVHRR measures a "brightness temperature" of 
the sea surface, and the actual temperature can be computed 
using physical radiation laws. The sea surface temperature 
is calculated from the brightness temperature by correcting 
for attenuation and emission of the atmosphere. The 
radiometer is calibrated on board the spacecraft by having 
the instrument view deep space (approximately 
2-degree-Kelvin black body), and the radiometer housing 
(monitored by thermistors) . The multi-channel simultaneous 
imaging of the sea surface forms the basis for the atmos- 
pheric corrections used for the AVHRR measurements. 

Deschamps and Phulpin ‘Ref. 41] proposed a 
mul tispect ral correction method which is the basis for the 
AVHRR measurement accuracy of sea surface temperature. The 
simultaneous solution of three different "window" 
measurements removes many non-linearities caused by 
atmospheric effects and approximations used zo calculate 
radiances. The cost of improving atmospheric corrections is 
an increase in the noise equivalent temperature of the 
system. 
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The AVHRR/1 (NGAA-6) uses two bands for analysis and 
atmospheric corrections. The analysis bands are channels 
three and four. Each channel has a different response to 

water vapor in the atmosphere. The variation of response of 
the two channels makes the corrections mors reliable. 
However, channel three operates very near the five 
micrometer wavelength where solar energy is about equal to 
the emitted energy of a 3 00-degree- Ke Ivin black body (sea 
surface). Therefore, channel three is a night channel. 
Channel four can be used day or night. 

Leirao et a 1 . [Ref. 42] discuss the use of satel- 

lite infrared imagery to locate the northern surface 
boundary of che Gulf Stream. They conclude that the maximum 
sea surface temperature gradient, is a good approximation to 
the location of the northern boundary. 

A problem with determining sea surface temperature 
gradients from space is the variability of the atmosphere. 
Maul [Ref. 43] discusses the effects of the atmosphere on 
ocean surface thermal gradients. In a moist atmosphere, a 
thermal gradient can be reduced to less than half of its 
actual value by the atmosphere. 

2. passive Micro wave R ad iometrv 

The latest instrument used for passive microwave 
observations of the earth is the Scanning Multi-frequency 
(channel) Microwave Radiometer (SMMR) . This instrument is 
operational on NIMBUS-7. Sea surface temperatures are 

calculated from brightness temperatures measured at 6.6 
gigahertz (4.54 cm wavelength). At this frequency, special 
resolution is generally considered to be 150 km. Test 
results of the NIMBUS-7 system indicated that resolution 
cells for the 6.6 GHz channel were 157x156 km [Ref. 44]. 
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the detrimental e f f sets 



Lipes [Ref. 45] referred to 
of rain and sunglint (reflection of the sun off of the water 
surface) on the SEASAT SMMR (a similar instrument is on 
NIMBUS-7). Rain and sunglint caused the SMMR derived sea 
surface temperature to be much higher than the actual sea 
surface temperature. Areas without sunglint or rain were 
more accurately measured by the SEASAT SMMR. 

In addition to problems with rain and sunglint, the 
SMMR is also susceptible to radio frequency interference 
(RFI) , and severe sidelobe contamination within 600 km of 
land [Ref. 46]. Improvements in the instrument design and 
in the ground-based processing design would reduce the side- 
lobe contamination problem for the SMMR [Ref. 47]. 

3. Satellite Radar Alt im eter 

Satellite altimetry presents an active approach to 
locating water mass boundaries and other density related 
anomalies in the oceans. The satellite measures height 
above the surface of the ear^h. These data, combined with 
the orbital parameters of the spacecraft, form the basis for 
measuring the height of the sea surface along the 
sub-satellite track. The radar altimeter maps the 

topography of the ocean along the satellite ground track, 
since the instrument is fixed at nadir and does not scan 
across track. The ocean topography is the difference 
between the sea surface and the geoii. Surface gecs J -rophic 
currents can be derived from the ocean topography generated 
by a radar altimeter. 

Altimeter height measurement uncertain ities exist 
due to atmospheric propagation and ionospheric effects, 
local deviations of the ocean surface topography from the 
geoid, and inaccuracies in the satellite orbit. Density 
gradients associated with geostrophic currents account for a 
small part of the topographic variations [Ref. 48]. 
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The SEASAT altimeter, accurate to 10 centimeters in 
absolute altitude measurement, was an excellent instrument 
for measuring the location of surface geostrophic currents. 
The G EOS -3 altimeter, accurate to 50 centimeters in absolute 
altitude, was capable of locating the slope associated with 
the Gulf Stream by using an averaging process [Ref. 49]. 

The altimeter is a promising instrument for the 
future. It is capable of detecting western boundary 
currents such as the Kurcshio, but the error of the height 
measurement is still larger than the signal generated by a 
weak current. A western boundary current will cause the sea 
surface to rise one meter over 100 kilometers. A 10 centi- 
meter accuracy of height measurement is required to measure 
such a slope of the sea surface accurately. When the slope 
of the sea surface is an order of magnitude less than 
Im/IOOkm, mere accurate height measurements are required. 

B. OTHER INSTRUMENTS FOR OBSERVING THE OCEANS 

This chapter has been directed toward those instruments 
which may be of importance tc the problem of predicting the 
mixed layer depth from sea surface temperature. The problem 
is much more complex than cause and effect between sea 
surface temperature and mixed layer depth. This report is 
limiting the scope of the investigation to ascertaining a 
relationship between sea surface temperature, thermocline 
gradient, and mixed layer depth, if indeed one exists. 

The problem can be expanded to include the effects of 
wind mixing and heat budget. The space technology needed to 
provide observations of wind speed, wave height, and temper- 
ature has been demonstrated on SEASAT. The SMMR can provide 
information on the wind speed, in addition to sea surface 
temperatures. Besides surface topography, the altimeter on 
SEASAT cculd provide information on significant wave height. 
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Additional instruments on SEASAT provided wind velocity 
(Scatterometer) , and surface wave direction and wavelength 
(SAR) . The Scatterometer (SCAT) and the Synthetic Aperature 
Radar (SAR) are active microwave imagers and not subject to 
many of the problems encountered by passive systems. 
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IV 



l^PERI MENTAL DESIGN 



A. OBJECTIVE 

The objective cf this investigation was no determine if 
a relationship between sea surface temperature, mixed layer 
depth, and thermocline gradient could be established. By 
reviewing the climatology of the area, and including the sea 
surface temperature patterns and sub-surface structure, an 
empirical approach to forecasting the extent of the mixed 
layer from surface observations and limited bathythermo- 
graphs was to be developed. 

Legeckis and Gordon [Ref. 50] found that 

sat el li 1 - e- derived sea surface temperature fronts 

corresponded to deep-reaching subsurface temperature 
structure. They also concluded that the observed mixed 
layer depth correlated with the surface temperature patterns 
observed for the Brazil Current and two warm core eddies. 

Blackstone and Whritner [Ref. 51 ] studied the persis- 
tence cf observed sea surface temperature patterns to a 
depth of 1000 feet off the coast cf southern California. 
Air-dropped Expendable Bathythermographs (AX3T) were 
employed to examine regions which displayed a strong sea 
surface temperature gradient in Defense Meteorological 
Satellite imagery. The AXBT's provided observations cf the 
vertical temperature structure associated with the 
observable sea surface temperature patterns. From analysis 
of temperature patterns at various depths, Whritner and 
^ Blackstone observed that the temperature gradient associated 
with sea surface temperature patterns was still recognizable 
at 1000 feet. 
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Leitao ef_ al. [Ref. 52] found that the 
surface temperature gradient (as observed from 
useful for approximating the northern surface 
the Gulf Stream. 



maximum sea 
space) was 
boundary of 



B. DATA ACQUISITION 

In November and December of 1980, the Naval Postgraduate 
School conducted the Acoustic Storm Transfer and Response 
Experiment (discussed in chapter 1) . As part of this exper- 
iment, Air-Dropped Expendable Bathythermographs ( AXBT) were 
placed along the same track on six days. There were three 
flights in November, covering the center track. In 
December, another three flights were conducted. The center 
track was again covered three times, but on the last two 
flights, tracks offset from the center also were flown. The 
offset tracks necessitated a reduction in the number of 
stations which could be taken along either the center track 
or the offset track [Ref. 53]. 

Satellite observations of sea surface temperature were 
collected from infrared imaging by the Advanced Very High 
Resolution Radiometer (AVHRR) on board the NOAA-6 
spacecraft. AXBT's were placed along a 1000-km crack (at 
55- km intervals). The data from the AXBT’s were digitized 
in five-meter increments of depth (see Figures 4.1 and 4.2), 
and the printed output was used to construct the vertical 
cross-sections of temperature used in this thesis. 

In addition to temperature soundings collected from the 
AXBT’s, soundings were made from a research vessel in the 
operating area using a Neil Brown conductivity temperature 
depth (CTD ) microprofiler [Ref. 56]. 
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Digi tiz ad 



Depths of 



Is otherms 



[Ref. 
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Figure 4.3 



FNOC Salinity Profile Along the Flight Track 
for 13 November. 
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c. 



DATA ANALYSIS 



1 . Data Plots 

As a means of beginning an analysis of the possible 
relationship between sea surface temperature patterns and 
corresponding patterns of mixed layer depth, vertical 
cross-sections of temperature for each flight track in 
ASTREX were constructed. In addition to these transects, 
sea surface temperature was plotted by position along the 
flight track as well as mixed layer depth. A correlation 
was noted visually between track position and both sea 
surface temperature and mixed layer depth. 

The measurements taken with the CTD were in the 
northwest corner of the ASTREX region. A Fleet Numerical 
Oceanography Center salinity profile for * he track (Figure 
4.3) shews very little variation of salinity profiles 
between stations 2 and 13. Since the variability of 

salinity in the Northeast Pacific is low, the soundings 
which were taken with the CTD were assumed to be reasonably 
representative of the area. 

2 • St a t ist ical P r oc ~d u re 

Sea surface temperatures and mixed layer depths were 
taken from Kilonski's [Ref. 57] digitized depths of the 
isotherms. The mixed layer depth at a location as reported 
by Kilcnski sometimes shewed a variation of one meter 
between the two data formats (digitized temperatures at 

five-meter increments vs. digitized depths of the 
isotherms) . 

To examine variations with position along the track, 
the outermost station on each flight was assigned as the 
origin. All other stations were listed by (nautical mile) 
distance from the origin. For example, station nine was the 



49 



outermost station on 19 November. It was assigned the 
origin position (0.0) for that flight track. The adjacent 
station (seventeen) is thirty nautical miles from station 
nine. Therefore, its position is listed as thirty (30). 

In addition to position, sea surface temperature, 
and mixed layer depth, the gradient at the thermocline was 
considered a necessary factor to include in the list of 
variables. Four depth intervals were selected for computing 
gradients in the thermocline. The variations of gradients 
with track position are plotted in Appendix A (November 
tracks) and Appendix B (December tracks) . 

From the vertical cross-sections, an initially small 
interval was considered necessary for computing the 
gradient. The smallest interval available from Kilonski's 
data was five meters. Since the mixed layer depths did not 
fall on multiples of five meters, the first level above the 
reported mixed layer depth was used as a reference level for 
computing temperature changes at the top of the thermocline. 
For example, if the mixed layer depth was 67 meters, the 
temperature at 65 meters depth was used as the temperature 
at the top of the thermocline. To smooth gradient varia- 
tions caused by this method of defining the temperature at 
the top of the thermocline, gradients were also computed 
over intervals of fifteen, fifty, and one-hundred meters. 

The computed temperature changes are listed in 
various tables associated with each flight track. The 
normalized values of gradient computed from each interval 
are listed in Appendix C. 

The data analysis was done using the interactive 
statistical computing system MINITAB [Ref. 58]. Basic 
statistics (means and standard deviations), correlation, and 
linear regression analysis were performed for each flight 
track. MINITAB correlation of data resulted in the 
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(r) as discussed 



computation of the correlation coefficient 
in Hiller and Freund [Ref. 59]. 

Three models of mixed layer depth were developed for 
each flight track. These models were simply the least 
squares regression cf mixed layer depth as a function of sea 
surface temperature, position, and gradient in the t her mo- 
di ne. If the correlation between sea surface temperature 
and mixed layer depth was weak, then the model of mixed 
layer depth as a function of sea surface temperature was 
ignored . 

Regression analysis in MINITAB yields a curve and a 
coefficient of determination as a measure cf the goodness of 
the curve-fit (the coefficient of determination is the 
square cf the correlation coefficient). To avoid confusion, 
each of the measures of correlation between various combina- 
tions of parameters will be identified whenever used. 

D. OBSERVATIONS BY FLIGHT 
1 . FI icht Track One 

The vertical temperature cross-section for flight 
track one is shown in Figure 4.4. Ten stations were used 
to construct this profile. Isotherms tended to be deeper 
nearer the coast than out at sea. 

A warm water pocket exists between stations four and 
two. This warm pocket persisted through the vertical extent 
of the soundings between these stations. 

Figure 4.5 shows both the eea surface temperature 
and mixed layer depth along the flight track. A possible 
correlation between sea surface temperature and mixed layer 
depth along the flight track is suggested by the variations 
of both with position along the track (r = 0.33) . 
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Station Number 




Figure 4.4 



Vertical Temperature Cross-section 
(1 5 November) . 



52 



Station Numbtr 



TEMPERATURE (*C> 




I 

Figure 4.5 Mixed Layer Depth and Sea Surface Temperature 
(15 November) . * 
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For flight one, a strong negative correlation (r = 
-0.818) exists for the gradient computed from the upper five 
meters of the thermocline. This suggests that the initial 
gradient at the thermocline weakens in the transition region 
(coastal waters) . The weakening of this initial gradient 
could also be due to the presence of a warm meander or eddy 
as suggested by the vertical temperature cross-3 action . 





Table IV 


lists, by st a 


t ion 


• 






1 . 


position re 


lative to the 


mcs 


t-sea ward buoy ; 






2. 


mixed layer 


depth; 










3. 


sea surface 


tern perat u re ; 










4 . 


temperature 


change in t 


he 


upper five meters 


of 


the 




thermocline 


(dT/5) ; 










5. 


temperature 


ch ang e i n 


the 


upper 15 meters 


0 f 


the 




t her mocline 


(dT/15) ; 










6 . 


temperature 


change in 


the 


upper 50 meters 


o f 


the 




thermocline 


(dT/50); and 










7. 


temperature 


change in 


the 


upper 100 meters 


of 


the 




thermocline 


(dT/100) . 










Table 


IV is the 


matrix which was manipulated to 


examine 


cor re 


lations between various 


elements in the problem. 


and 


was 


a source of 


values in r 


egression analysis fo 


r* 


model 



equat ions. 

Attempts to fit several equations for mixed layer 
depth were made using various measurable quantities. The 
results are listed in Appendix D. Table V lists "-he model 
variables used and provides the computed coefficient of 
determination for each model. The most striking relation- 
ship occurred for a fit of mixed layer depth as a function 
of track position, sea surface temperature, and gradient 

from 100 meters. The coefficient of determination for the 
fit was 0.88. The correlation (r) between MLD and the func- 
tion of the three variables was 0.94. Position and 
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T ABLE IV 



MLD, SST and Gradients at the Thermocline, 15 November 1980. 



Buoy 


Position 


MLD 


SST 


dT/5 m 


dT/ 15m 


dT/50m 


dT/IOOm 


2 


54 0 


32 


13.81 


1 . 15 


2. 86 


3.90 


4.65 


3 


480 


69 


1 5.88 


1.41 


2.69 


4.62 


5.92 


4 


420 


59 


1 4.86 


1. 13 


3. 14 


5.35 


6.22 


5 


36 0 


45 


14.43 


1. 56 


2. 58 


4.01 


5.06 


6 


30 0 


62 


1 4.32 


1.70 


3.64 


5.27 


6.49 


7 


24 0 


51 


1 3.4 1 


2.52 


3. 36 


4.22 


5.52 


8 


180 


6 1 


1 3.27 


1.94 


3. 33 


4.49 


5.66 


9 


120 


51 


1 2.52 


2. 48 


3. 04 


4.01 


5.30 


1 0 


60 


53 


1 2. 33 


2.67 


3.66 


4.74 


5.49 


1 1 


0 


54 


1 2.09 


1.99 


2. 88 


4.09 


4 . 79 



Not 9S I 

1. ~Bucy column indicates station number in vertical 

temperature cross-sections. 

2. Position column lists the range (in nautical miles) 
from the most seaward buoy. 

3. dT/5m is the temperature change over the first five 
meter interval in the thermccline; dT/15 is the 
temperature change over the first fifteen meter 
interval in the thermocline; dT/50 and dT/100 
follow in the same manner. 



TABLE V 



Models of ML D for 15 November. 



Model Variables 


Coefficient of 


SST 


0.111 


SST, Position 


0.818 


SST, Position, 
Gradient (100) 


0.880 



on 



Notes : 

1. SST is sea surface temperature. 

2. Position refers to station position 
relative to the most seaward station. 

3. Gradient (100) refers to the gradient 
computed from the upper 100 meters of 
the thermccline. 

4. The coefficient of determination will equal 
one for a perfect fit of the curve to data. 
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temperature are easily measured by satellite. Gradient then 
becomes the primary factor of uncsrt ainit y . 

2 - Fl ioh t Track Tw o 

The vertical temperature cross-section for 17 
November is shown in Figure 4.6. Twenty-three stations were 
used to construct this transect. The profile is interesting 
for three reasons: 

1. This flight has no observational "holes 1 ' along the 
track ; 

2. Satellite imagery from NOAA-6 is available for a 
portion of the track; 

3. The vertical temperature structure of a 
satellite-observed sea-surface- temperature anomaly can 
be studied. 

At the top of the thermocline, the gradient appears 
strong between stations 13 and 16 (isotherms are densely 
packed). Beginning at station nine, the isotherms begin to 
diverge. This indicates a weakening of the gradient in the 
upper five to fifteen meters of the thermocline. As with 
flight one, the initial gradient (from five meters) is 
inversely correlated with position along this track. The 

initial strength of the gradient at the top of the thermo- 
cline decreases towards the coast. 

The location of the divergence of the isotherms at 
the top of the thermocline (station 9) corresponds to the 
location of the change from the Central Subarctic Domain to 
the Transitional Domain as presented by Dodimead ei_ al . 

Figure 4.7 is a plot of sea surface temperature and 
mixed-layer-depth variation with position. From station 
nine to station twenty-one, the sea surface temperature 
increases almost linearly with position, while mixed layer 
depth fluctuates around sixty meters. Satellite imagery 
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Station Numbtr 




Figure 4.6 



Vertical Temperature Cross-section 
(17 November^ . 
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identified a warm meander between station four and station 
twenty-two (Figures 4.8 and 4.9). This shows up as a rapid 
increase in sea surface temperature aver a (longer) distance 
from station four to station three. The sea surface 
temperature reaches a maximum in the vicinity of station 
three, and then falls rapidly across the southeast wall of 
the meander. The mixed layer depth also reaches a maximum 
at station three, and then falls rapidly across the 
southeast wall. This correlation between sea surface 
temperature patterns and vertical structure in the vicinity 
of frontal positions was also noted by Legeckis and Gordon 
(Ref. 33), cited earlier. 

Examinations of this track based on domains and 
satellite imagery requires that the track be split into 
three separate sections. The first section, in the Central 
Subarctic Domain, includes stations thirteen through 
sixteen. The Transitional Domain is then split using the 
satellite imagery. One section extends from station nine to 
station twenty-one. The warm meander then constitutes the 
final section and includes stations four to twenty-four. 

Table VI provides position, mixed layer depth, sea 
surface temperature, and gradients in the the r mociine for 
each station occupied for flight track :wo. As with flight 
one, this table constitutes the matrix used to examine 
relationships between various elements which may lead to a 
prediction of mixed layer depth. 

Correlations of sea surface temperature and mixed 
layer depth along the entire track yield low values for the 
correlation coefficient (r) . For the entire track, there is 
little correlation between sea surface temperature and mixed 
layer depth (r = 0.17). Mixed layer depth as a function of 
position, sea surface temperature, and gradient from the 100 
meter ther mociine is better (r = 0.84). However, along 
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Station Numbir 



TEMPERATURE t*C> 




Figure 4.7 Mixed Layer Depth and Sea Surface Temperature 
(17 November) . 
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Figure 4.9 



Flight Track. Buoy Numbers and Meander Location 
in Figure 4.3. 
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T ABLE VI 



MLD , SST and Gradients at the Thermocline, 17 November 1980. 



Buoy 


Position 


MLD 


SST 


dT/5 m 


dT/15m 


aT/50m 


dT/IOOm 


2 


63 Q 


45 


1 3.76 


1.33 


2.35 


3.40 


3.97 


3 


570 


70 


1 5.80 


1. 51 


2. 82 


4.77 


b. 53 


4 


510 


67 


1 4. 83 


1.48 


3. 71 


5. 01 


6.04 


5 


45 0 


54 


14.46 


1.51 


3. 47 


4.95 


6.11 


6 


39 0 


61 


1 3.97 


1. 64 


2. 93 


4.09 


5.33 


7 


330 


54 


1 3.06 


1.69 


3. 25 


4.33 


5.70 


8 


270 


70 


1 2.87 


2.65 


3. 35 


4.64 


5.30 


9 


210 


64 


1 2. 14 


1. 19 


2. 96 


4.36 


4 .95 


1 0 


150 


58 


11.95 


1. 67 


3. 70 


4.76 


5.52 


1 1 


90 


58 


1 1.63 


2. 58 


3. 39 


4.73 


5 . 19 


1 2 


30 


58 


10.61 


1.85 


2. 63 


3.87 


4.23 


13 


0 


56 


10.13 


2.47 


3. 09 


3.90 


4.25 


14 


60 


48 


10.72 


1.71 


2. 68 


3.62 


4.24 


1 5 


120 


56 


11.71 


2. 58 


3. 36 


4.47 


5.01 


1 6 


180 


58 


1 1.83 


2. 08 


3.72 


4.82 


5.25 


17 


24 0 


53 


1 2.65 


1.75 


2.38 


3.90 


4 .95 


1 8 


300 


59 


1 2.90 


1.43 


2. 88 


4.07 


5.64 


1 9 


360 


5 3 


13.86 


2.66 


3. 63 


5.06 


6.22 


20 


420 


54 


1 4.43 


1 . 16 


3.55 


5.20 


6 .24 


2 1 


48 0 


58 


14.5 1 


1.83 


3. 84 


4. 84 


6.01 


22 


540 


53 


1 5.24 


1.80 


3. 50 


5.70 


6.70 


23 


600 


69 


15.59 


0.93 


2. 64 


4.66 


6 . 30 


24 


66 0 


24 


1 3.54 


1. 08 


1.9 5 


4.28 


5.35 



shorter track segments, and particularly across the meander, 
mixed layer depth and sea surface temperature show more 
correlation than for the entire -crack. Across the meander, 
the correlation coefficient for mixed layer depth vs. sea 
surface temperature is 0.86. This was not the case in 
either the Central Subarctic Domain (r = 0.40) or in the 
other part of the Transitional Domain (r = -0.53). 

The mixed layer depth and sea surface temperature 
did not correlate to any siginificant degree on 17 November. 
Therefore, only two models are proposed for this flight 
track. As stated previously, breaking the track into 
segments yields a much higher correlation between otherwise 
uncorrelated parameters. Table VII presents the models for 
the entire track only. 
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TABLE VII 



Models of MLD for 17 November. 



Model Variables Coefficient of Determination 



SST , Position 

SST , Position, 
Gradient (15) 



SST . 
Grad 



Position, 
ient (100) 



0.615 

0.627 

0.707 



Not e ! 

1. Two gradient models were computed 
because the gradient from 15 mere 
appeared to be more highly corral 
wiuh track position than the grad 
computed from the upper 100 meter 
the' thermocline. 



r s 

a ted 

TIT 

s of 



The higher correlation between sea surface tempera- 
ture and mixed layer depth in the warm meander may provide a 
means of estimating the depth of the mixed layer in such 
anomalies through observations of sea surface temperature. 

During the interval between flight one and flight 
two, one storm had traversed the area and another had begun 
to approach the most western stations on the track. A 
slight deepening of the 10-degree through 11-degree 
isotherms occurred. The average mixed layer depth changed 
from 53.7 meters on 15 November to 57.0 meters on 17 
Nov em ber . 



3 . Flight Track Three 

Figure 4. 10 is the 
with flight track three, 
construct this transect, 
the data along the track, 
similar to flight track two 



vertical cross-section associated 
Fourteen stations were used to 
Although there are some holes in 
this profile presents a picture 
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The gradient at the top of the thermocline diverges 
beginning at station nine, just as in the previous secticns. 
Again, the gradient computed from the first five meters of 
the thermocline was inversely correlated with position (r - 
-0. 766) . The meander discussed for flight two still can be 
located in this profile between stations four and two. 

The sea surface temperature and mixed layer depth 
variations with position are plotted in Figure '4.11. Table 
VIII lists the values for mixed layer depth, sea surface 
temperature, and gradients alone the flight track. Between 
station nine and station six (Transitional Domain) both sea 
surface temperature and mixed layer depth vary almost 
linearly. Over the entire track length, however, the 
correlation between the two is negligible (r = 0.10) . From 
station twenty-one to station one (across the meander) the 
correlation between mixed layer depth and s-i-e surface 
temperature is much better (r = 0.51). 



A 


front 


move d 


into the 


flight 


track area on 


Nov em ber , 


and 


beca me 


st atior.a 


ry over 


the area on 


November . 


The 


stati cna 


ry front 


caused 


poor conditions 



satellite observations due to clouds. Small changes in sea 
surface temperature par ter ns could be expected from -he 
front, but the basic structure remained unchanged from 17 
November. The average mixed layer dsp^h for 19 November 
remained at 57.0 meters. The position and structure of ^ he 
meander along the track remained unchanged except for the 
development of a peak temperature exceeding 16 degrees 
Celsius at station three. 

Table IX summarizes the models examined for flight 
track three. The correlation coefficient computed for sea 
surface temperature versus mixed layer depth was extremely 
low (r = 0.01) . Therefore, the model equation for mixed 
layer depth as a function of sea surface temperature was 
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Station Numbtr 




Figure 4.10 Vertical Temperature Cross-section 
(19 Noveaber) . 
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Figure 4.11 Mired Layer Depth and Sea Surface Temperature 
(19 November) . 
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TABLE VIII 



MLD, SST and Gradients at the T hennas line , 19 November 1980. 



Buoy 


Position 


MLD 


SST 


dT/5 m 


dT/ 15m 


dT/50m 


dT/IOOm 


1 


48 0 


30 


13.17 


1 . 15 


2.86 


3.90 


4.65 


3 


36 0 


67 


16.04 


1.41 


2. 69 


4.62 


5 .92 


4 


300 


58 


14.8 1 


1. 13 


3. 14 


5.35 


6 . 22 


5 


24 0 


59 


14.43 


1.56 


2. 58 


4.01 


5.06 


6 


180 


54 


1 3.95 


1.70 


3.64 


5.27 


6.49 


7 


120 


53 


13.14 


2.52 


3.36 


4.22 


5.52 


8 


60 


62 


1 3.00 


1.94 


3. 33 


4.49 


5.66 


9 


0 


66 


1 1.98 


2.48 


3. 04 


4.01 


5.30 


17 


30 


64 


1 2.79 


2.67 


3. 66 


4.74 


5.49 


18 


90 


61 


1 2.95 


1.99 


2.88 


4.09 


4 .79 


20 


210 


64 


14.32 


1. 15 


2. 86 


3.90 


4.65 


2 1 


27 0 


54 


14.35 


1.41 


2. 69 


4. 62 


5 .92 


22 


330 


43 


1 5.80 


1 . 13 


3. 14 


5.35 


6.22 


23 


390 


59 


15.48 


1.56 


2.58 


4.01 


5.06 


24 


45 0 


54 


13.51 


1.70 


3. 64 


5. 27 


6.49 



omitted. Instead, a second gradient model was introduced as 
with flight track two. The gradient from 100 meters was 
continued as an element in a separate model because this 
gradient was either the most influential or second most 
influential of all gradients used in model equations. 

The best fit of any model for mixed layer depth was 
the model based only on sea surface temperature and track 
position. Both surface temperature and track position are 
variables which can be derived from a satellite image of the 
sea surface. 

4 . FI ight Track Four 

The vertical cross-section of temperature for flight 
track four (Figure 4.12) was constructed from 21 stations. 
It is the most complete transect except for 17 November. A 
definite cooling has occurred between 17 November and this 1 
December track. Surface temperatures dropped by about one 
degree Celsius (the average sea surface temperature dropped 
from 13.1 degrees Celsius on 17 November to 12.0 degrees 
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TABLE IX 



Models of MLD for 19 November. 



Model Variables Coefficient of Determination 



SST , Position 


0.615 


SSI* Position, 
Gradient (5) 


0.547 


SST, Position, 
Gradient (100) 


0.517 



Not e : 

1. Two gradient models were computed 
because the gradient from 5 meters 
appeared to Be more highly correlated 
with track position than the gradient 
computed from the upper 100 meters of 
the thermocline. 



Celsius on 1 December). The gradient in the upper five to 
fifteen meters of the thermocline is not as strong as along 
the November flight tracks. Although still negatively 
correlated with track position, the correlation between 
gradient and track position is much weaker for this track 
{r = -0.16) than for similar gradients computed from the 
upper five meters of the thermocline. The meander is still 
recognizable, but its vertical profile is much less 
pronounced than during the November flights. 

The gradient computed from the upoer five meters of 
the thermocline showed the strongest correlation with posi- 
tion for all gradients computed along flight track four. As 
noted previously, the correlation between the five meter 
gradient and position has steadily decreased since the 
first flight. However, the five meter gradient is the most 
highly correlated of any gradient with position on each of 
the four flights. 

Mixed layer depth doesn’t show much relationship to 
track position. Sea surface temperature dees increase along 
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the track, but the change across the meander is not as 
significant as in past flights. Figure 4.13 displays the 
variation with position for both mixed layer depth and sea 
surface temperature. 



TABLE X 

MID, SST and Gradients at the Thermocline, 1 December 1980. 



buoy 


position 


MLD 


SST 




dT/5 m 


dT/15m 


dT/50m 




dT/ 100m 


1 


690 


51 


1 2.98 




1. 15 


2.86 


3.90 




4.65 


4 


510 


61 


14.46 




1.41 


2.69 


4.62 




5.92 


5 


45 0 


64 


1 3. 30 




1. 13 


3. 14 


5.35 




6.22 


6 


39 0 


69 


1 2.92 




1.56 


2. 58 


4.01 




5.06 


7 


330 


72 


1 1.74 




1.70 


3. 64 


5.27 




6.49 


8 


27 0 


72 


1 1.82 




1. 94 


3.33 


4.49 




5.66 


9 


21 0 


61 


11.12 




2.52 


3. 36 


4.22 




5.52 


10 


150 


72 


1 0.72 




1.94 


3. 33 


4.49 




5.66 


1 1 


90 


78 


1 0.77 




2.52 


3. 36 


4.22 




5.52 


. 12 


30 


70 


9.13 




2. 48 


3.04 


4.01 




5.30 


1 3 


0 


74 


8.94 




2.67 


3.66 


4.74 




5.49 


14 


60 


72 


9.5 1 




1.99 


2. 88 


4.09 




4.79 


1 5 


120 


69 


1 0.32 




1. 15 


2. 36 


3.90 




4.65 


17 


24 0 


70 


1 1.42 




1.41 


2.69 


4.62 




5 .92 


1 8 


300 


78 


1 1.70 




1 . 13 


3. 14 


5.35 




6 . 22 


19 


36 0 


74 


1 2.98 




1. 56 


2.58 


4.01 




5.06 


20 


420 


80 


1 3.25 




1.70 


3.64 


5.27 




6.49 


2 1 


48 0 


69 


13.54 




1.99 


2. 38 


4.09 




4.79 


22 


540 


66 


1 4.24 




2. 67 


3. 66 


4.74 




5.49 


23 


600 


59 


14.27 




2.48 


3. 04 


4.01 




5.30 


24 


66 0 


50 


1 3.06 




2.67 


3. 66 


4.74 




5.49 




Table X 


li 


st s the 




factors 


which 


were e 


xamined for 


flight 


four. 


This 


f light 




was unusual in 


that la 


r 


ge changes 


from p 


revicus fligh 


ts occu 


■r 


red in 


the corr 


elation 


s 


• 




Both sea surface 


im 


empe rat 


urs and 


mix ed 


l 


ayer depth 


had higher corr 


elations wi 


r 


h track 


position than 


on previous 


flights. The 


CO 


rr elati 


on of 


graaien 


ts wit 


h 


posit ion- 


declined in all 


four cases 




examine 


d. Mix 


ed layer 


depth was 


mor e h 


ighly cor 


related wit 


h 


sea surface t 


emperat 


u 


re than on 



previous flights. 
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Station Humber 



f 



r 



T 




Figure 4.12 



Vertical Temperature Cross-section 
(1 December) . 
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Station Numbtr 



i 



TEMPERATURE CC> 




Figure 4.13 Sized Layer Depth and Sea Surface Temperature 
(1 DeceniDer) . 
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Three models are presented for flight four and 
summarized in Table XI. Since the correlation of sea 
surface temperature and mixed layer depth was the highest 
observed of all flights, i t is included as a model. Mixed 
layer depth as a function of track position, sea surface 
temperature and thermocline gradient showed the best 
correlation when the gradient computed from 100 meters was 
used. 



TABLE XI 

Models of ML D for 1 December. 



Model Variables 



Coefficient of Determination 



SST 



0.175 



SST, Position 



SST. Position 
Gradient (100 



f 



0.609 

0.663 



Not e : 

1. One gradient model was computed 

because the Gradient from 100 meters 
appeared to Be more highly correlated 
with track position than any other 
gradient . 



5. Flight. Track Five (Center and South) 

The fifth flight added another dimension to x he AXBT 
analysis. The southern track, offset by 60 nautical miles 
from the center track, provided additional information about 
the extent of the warm meander discussed previously. Figure 
4.14, constructed from eleven stations, is the vertical 
Temperature Cross-section for the center track on 3 
December. Figure 4.15, constructed from eight stations. 
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shows the profile for the southern track. Table XII lists 
the values for mixed layer depth, sea surface temperature, 
and gradients along the center flight track. Table XIII 
lists similar data for the southern track. 

The vertical temperature profiles for the two tracks 
are similar in some respects. The meander is definitely 
present on the center track. The southern track displays 
some similar subsurface structure in the area where the 
meander should be located (between stations forty-one and 
forty-three). Although station forty-two should be inside 
the warm meander, the surface temperature at the station is 
13.5 degrees Celsius, about one-half degree too low to be a 
part of the meander. The closed warm pocket of water is not 
present in the southern track (the center track has a closed 
pocket of 14 degree Celsius water). 

The meander/eddy influence appears to be present to 
about 290 meters in both tracks. The southern track does 
not have the pronounced structure shown along the center 
track. Although some signs of the meander remain in the 
southern track, the meander appears to have formed a closed 
cir cu lat ion. 

Along the center flight track, sea surface 
temperature and mixed layer depth (Figure 4.16) are not 
highly correlated (r = 0.38). Across the eddy, however, th a 
correlation coefficient for sea surface temperature versus 
mixed layer depth increases to 0.88. The southern track 
(Figure 4. 17) also shows an excellent correlation between 
sea surface temperature and mixed layer depth ( r = 0.84). 
The corrlea tion between these two variables is less (r = 
0.66) across the track segment which coincides with the 
eddy. It appears that the meander has formed into an eddy 
which dees not extend to the southern track. 
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Station Number 




Figure 4.14 Vertical Temperature Cross-section 
(3 Dec, Center) . 
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Station Number 




Figure 4.15 Vertical Temperature Cross-section 
(3 Dec, South ) . 
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Figure 4.16 



Mixed Layer Depth and Sea Surface Temperature 
(3 Dec, Center). 
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Figure 4.17 Mixed Layer Depth and Sea Surface Temperature 
(3 Dec, South ) . 
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Most of the southern track is in the Transitional 
Domain. The Central Subarctic Domain has influence up to 
station thirty-five, after which the track is under the 
influence of the Transitional Domain. Changes apppear to be 
linear along the track within the Transitional Domain. 

Winds increased from 11 meters/sec (24 lets) to 14-16 
meters/sec (30-35 kts) on 1 December. These wind speeds 
were sustained through 2 December and ■‘■hrough most of 3 
December [Ref. 61 ]. These cold northwesterly winds lowered 
the sea surface temperature along the track, but caused no 
major changes to the vertical Temperature Cross-section. 

The average mixed layer depth for the center track 
was 66.4 meters on 3 December. The southern track averaged 
64.6 meters. There was a decrease in the average mixed 
layer depth along the track since the previous flight on 1 
December. In comparing the two days, a large decrease in 
mixed layer occurred at station five (21 meter change). At 
first glance this appears to be an instrument error. 
However, the mixed layer depth profile for flight track one 
(Figure 4.5) is very similar to the profile for flight track 
five. On both days station five had a shallow mixed layer 
which deepened considerably across the eddy. 

When overlaid, the mixed layer depth profiles and 
the sea surface temperature profiles for tracks one and five 
are remarkably similar. The mixed layer depth for flight 
five is deeper than for flight one, but the variations with 
position are very similar. 3etween. stations ten and two, 
the sea surface temperature is highly correlated (r = 0.98) 
with position, and the mixed layer dep^h shews some 
correlation (r = 0.54). 

Four model equations were used for each track on 
flight five. The models are summarized in Table XIV for the 
center track, and in Table XV for the southern track. On 
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TABLE XII 



MLD , SST and Gradients at the Thermosline, 3 December 1980. 



Buoy 


Position 


MLD 


SST 


dT/5m 


dT/1 5m 


dT/50ra 


d?/1 00m 


46 m 


690 


45 


1 2.47 


0. 04 


1. 28 


3. 58 


4.48 


3 


570 


70 


1 4.64 


1 . 05 


2. 39 


4. 67 


5.97 


4 


510 


62 


1 4. 16 


1.35 


2. 73 


4.87 


5.46 


5 


450 


42 


1 3. 38 


1.41 


2. 17 


3.35 


4 .51 


6 


390 


72 


1 3.06 


1. 42 


2. 44 


4.20 


5. 20 


7 


33 0 


69 


1 1.77 


0. 81 


2. 07 


3.01 


4.14 


8 


27 0 


83 


11.42 


1.53 


2. 66 


3.71 


4. 36 


9 


21 0 


77 


1 0.99 


1.80 


2. 26 


3.24 


3.37 


10 


150 


62 


10.02 


2. 28 


2.71 


3.30 


3.82 


12 


30 


66 


8.86 


0.78 


1 . 43 


2. 24 


2.43 


1 3 


0 


82 


8.46 


1.42 


1.38 


2.20 


2.49 



TABLE XIII 



MLD, SST and Gradients at the Thermocline, 3 December 1980 

(South) 



Buoy Positi 


on MLD 


SST 


dT/5 m 


dT/ 15m 


dT/50m 


dT/IOOm 


34 


0 


80 


9.40 


1. 23 


2. 56 


3.23 


3.42 


35 


60 


86 


9.53 


1. 43 


1 . 70 


2. 40 


2 . 62 


38 


24 0 


70 


11.44 


1.69 


1. 38 


3. 08 


3.55 


39 


30 0 


64 


1 2.04 


1.00 


2. 53 


3.34 


4.51 


4 1 


420 


74 


1 3.27 


0.92 


2. 60 


3.43 


4 . 89 


42 


48 0 


43 


1 3.46 


1.79 


3. 61 


4.20 


4.76 


43 


54 0 


46 


1 3.89 


1.62 


2. 56 


4.20 


5.12 


44 


600 


54 


13.65 


1. 07 


2. 10 


4.09 


5.22 


.f 

(D 


south t 


rack. 


the gradient 


compute 


d from 


100 meters 


produced the 


worst 


fit o 


f any 


model of 


mixed 


layer as a 


func 


tier, of 


sea surface 


t empe r 


at ure , 


track position and 


ther 


mocline 


gradien 


t. The 


gradi 


ent from 


100 meters produced 


(i) 

,r: 


best fit 


for t 


he cent 


er track model 


. For 


comparison. 


Ct- 

tr 

(D 


s cut hern 


track 


modal 


u sing 


the grad 


ient computed from 



100 meters is included. 
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TABLE XIV 



Models of MID for 3 December, Center Track. 

Coefficient of Determination 
0.146 
0.404 

0.44 1 

0.523 

TABLE XV 

Models of MLD for 3 December, South Track. 

Coefficient of Determination 
0.703 
0.716 

0.86 1 

0.716 

6. Flight Track Six (Center and North) 

Flight six examined both the center track (Figure 
4.18) and a parallel track 60 nautical miles to the north 
(Figure 4.19) with eleven stations each. The values of 
mixed layer depth, sea surface temperature, and gradient 
associated with each station are given in Table XVI for the 
center track, and Table XVII for the northern track. The 
temperature structure along the center track resembles 



Model Variables 
SST 

SST , Position 

SST. Position, 
Gradient (50) 

SST, Position. 
Gradient (100) 



Model Variables 
SST 

SST, Position 

SST, Position, 
Gradient (15) 

SST, Position. 
Gradient (100) 
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previous flight tracks. The isotherms begin no diverge 
after station nine. The gradient computed from the upper 
five meters of the thermocline is negatively correlated with 
position (r = -0.27). 

The eddy seems to have diminished both in length 
along the center track and in depth. The northern track 
does not appear to traverse the eddy. 

The sea surface temperature and mixed layer depth 
variation with position along the center track (Figure 4.20) 
show the best correlation in any flight along the center 
track. Both sea surface temperature and mixed layer depth 
varied linearly with position up to station three. Sea 
surface temperature correlation wi'-h position along the 
entire track was better (r = 0.92) than along the shorter 
track across the eddy (r = -0.80) , cr along any previous 
center flight -crack. Mixed layer depth was mors highly 
correlated with position across the eddy (r = -0.95) than 
along the entire track (r = -0.83). 

Across the eddy on the center track, sea surface 
temperature and mixed layer depth again are highly 
correlated (r = 0.90). Over the entire track length, the 
relationship between sea surface temperature and mixed layer 
depth is nor as strong (r = -0.55). Of all flights, this 
was the best overall correlation of sea surface temperature 
and mixed layer depth for the center track. 

The northern track (Figure 4.21) shows a 
relationship similar to that along the center track between 
sea surface temperature and mixed layer depth. The eddy is 



not apparent in 


the sea 


surface temper 


ature profile, but 


does present a 


f am ili ar 


profile of 


mixed layer depth. 


Across the eddy 


position 


on the north 


track, mixed layer 


depth is highly 


correlated 


with position 


(r = - 0 . 97) . This 


same relationsh 


ip holds 


true for fc he 


center track (r = 



-0. 95) . 
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Station Number 




Figure 4.18 Vertical Temperature Cross-section 
(5 Dec, Center) . 
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Station Numbtr 




Figure 4.19 



Tertical Temperature Cr oss-section 
(5 Dec r North) . 



33 



Station Number 



temperature (*C) 




Figure 4.20 



aired Layer Deprh and Sea Surface Temperature 
(5 Dec, Center) . 
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Station Number 



temperature cc> 
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Figure 4.21 Sized Layer Depth and Sea Surface Tempera 
(5 Dec, North ) . 
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A subsurface feature developed over the course of 
the experiment, and became very obvious in the northern 
track on 5 December. At station 58, the isotherms reach a 



maximum depth. 


A 


co rrespo ndin g 


excursion 


to depth for 


the 


same isotherms 


(6 


and 7 degree 


Celsius) 


occurs along 


^he 


center track at 


station nine. 


In looking 


back through 


the 



previous flight cross-sections, this same excursion is found 
to occur on flights four and five. 

Lundell [Ref. 62] reported a submerged sound channel 
in the vicinity of station nine. The Kiicnski digitized 
bathythermographs indicated a slight inversion in this area 
(of the order 0.3 degrees Celsius). This unusual structure 
occurs at the boundary between two Domains: the Central 
Subarctic and Transitional. Along the center track, the 
maximum depth excursion of the six and seven degree 
isotherms shifted to the west between flight two and flight 
six . 

This submerged feature is not driven by surface 
conditions as in occurs well below both thermoclir.e and 
halocline. This appears to be some type of shift in lower 
domain boundaries. A time series over one or twc months, 
between stations seven and ten, might have provided some 
insight into the lower layer variability at this domain 
boundary . 

As with the southern track on flight five, the 
northern track data on flight six showed the worst fit by a 
model using the gradient computed from 100 meters. The 
gradient from fifty meters produced the best fit on this 
track. The gradient from five meters produced the best fit 
on the center track. The model parameters and coefficients 
of determination are presented in Tables XVIII and XIX. 



86 



TABLE XVI 



MLD , SST and Gradients at the Thermocline, 5 December 1980. 



Buoy 


Position 


MLD 


SST 


dT/5m 


d T / 1 5 m 


dT/50m 


dT/IOOm 


1 


690 


42 


1 2.63 


0.84 


1.17 


3. 19 


4.21 


2 


63 0 


58 


13.11 


1.77 


3. 17 


3.79 


4.40 


3 


570 


70 


14.40 


0.89 


2. 10 


4.25 


5.83 


4 


510 


72 


13.8 1 


1 . 79 


3. 05 


4.60 


5. 44 


6 


390 


78 


1 2.63 


1.77 


3. 05 


4.00 


4 . 33 


7 


330 


75 


1 1.77 


1. 10 


2. 18 


3.40 


4 . 30 


8 


27 0 


78 


1 1.23 


1. 39 


2. 45 


3. 67 


4.25 


9 


210 


75 


1 0.88 


1. 11 


2. 46 


3. 42 


3.42 


1 0 


150 


78 


9.94 


1. 76 


2. 57 


3.47 


3.75 


12 


30 


82 


8.97 


1.60 


2. 26 


2.91 


2.33 


13 


0 


83 


8.43 


1 . 46 


1. 83 


2.02 


2.32 



Note: 

1. A slight inversion existed at station 12 between 
135 meters and 20 0 meters. 



TABLE XVII 



ML D , SS? and Gradients at the Thermocline, 5 December 1980 

( North) 



Buoy 


Position 


MLD 


SST 


dT/5 m 


dT/ 15m 


dT/50m 


dT/1 00m 


54 


0 


70 


8.8 1 


1 . 22 


1.59 


2. 14 


2.37 


55 


60 


83 


9.94 


1 . 36 


2. 32 


3.39 


3.69 


56 


120 


74 


10. 15 


2.03 


2. 68 


3. 65 


3.77 


57 


18 0 


77 


1 0.37 


2.24 


2. 89 


3. 79 


4.23 


58 


240 


6 4 


1 0. 96 


0.91 


2. 1 6 


2. 39 


2.69 


59 


300 


59 


1 1.20 


1 . 70 


2. 3 1 


3.28 


4.49 


60 


36 0 


70 


1 1.93 


0. 95 


1 .97 


3.03 


3.52 


6 1 


420 


69 


1 2.60 


0.92 


2. 28 


3.43 


4 . 39 


62 


480 


69 


1 2.73 


1.64 


3. 42 


3. 96 


4. 62 


64 


600 


50 


1 2.84 


1. 86 


3.01 


3.87 


4.71 


65 


66 0 


45 


1 2.63 


1. 12 


1 . 44 


2. 28 


3.59 



E. DYNAMIC HEIGHT COMPARISONS, 17 NOVEMBER 

Emery and O'Brien [Ref. 63] proposed a method of infer- 
ring salinity from temperature or depth in the North Pacific 
Ocean. The method relies on mean salinity values computed 
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TABLE XVIII 



Models of 


MLD for 5 


December, Center Track. 


Model Variables 


Coeff ici 


ent of Determination 


SSI 




0.298 


SST , Position 




0.962 


SST . Position, 
Gradient (5) 




0.977 


SST, Position, 
Gradient (100) 




0.963 



Models 


TABLE XIX 

of MLD for 5 December, North Track. 


Model Variables 


Coefficient of Determination 


SST 


0.370 


SST, Position 


0.846 


SST, Position, 
Gradient (50) 


0.872 


SST. Position, 
Gradient (100) 


0.854 



by five degree squares. In the test region, the mean 
salinity curves are nearly constant. Bernstein et al . 
[Ref. 64] used this method *:o estimate dynamic heights in 
the Western North Pacific with good success. It was decided 
to employ both a single recent salinity sounding in the area 
(from a CTD cast) , and a Fleet Numerical Oceanographic 
Center salinity climatology for November, as a basis for 
calculating and comparing dynamic heights for flight -rack 
two (17 November). 
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TABLE XX 



Dynamic Heights Relative to 350 meters (17 Nov). 



BODY 


POSITION 


CONST-S 


FNOC-S 


13 


0 . 


0.618 


0.629 


12 


30. 


0.643 


0.6*0 


14 


60. 


0.638 


0. 644 


1 1 


90. 


0.654 


0. 658 


15 


120. 


0.654 


0.658 


10 


150. 


0.656 


0.657 


16 


180. 


0.665 


0.663 


9 


210. 


0.695 


0.690 


17 


240. 


0.705 


0.698 


8 


270. 


0.713 


0.698 


18 


300. 


0.717 


0.690 


7 


330. 


0.714 


0.631 


19 


360. 


0.742 


0.699 


6 


390. 


0.741 


0.76 1 


20 


420. 


0.743 


0.753 


5 


450. 


0.759 


0. 768 


21 


480. 


0.766 


0.774 


4 


510. 


0.777 


0. 779 


22 


540. 


0.763 


0. 747 


3 


570. 


0.321 


0. 32 1 


23 


600. 


0.817 


0.316 


2 


630. 


0.737 


0.735 


24 


660. 


0.704 


0. 697 



Not es : 

1. Heights in dynamic meters. 

2. Const- S refers to sinale salinity profile for entire area. 

3. FNOC-S refers to November salir.it y‘ climatology from the 
Fleet Numerical Oceanographic Center. 



Consideration then was given to plotting the dynamic 
topography along the track. Since only one salinity profile 
was used to calculate the dynamic height at all stations, 
the resulting variations in dynamic heights along the track 
were due to temperature variations only. The reference 
level was selected as 350 meters because this depth repre- 
sented the useful limit of temperature data available from 
the digitized report. 

As a means of comparing dynamic height variation with 
position and other factors of interest such as temperature 
gradient at the thermocline or mixed layer depth, the Fleet 
Numerical Oceanographic Center climatology for the flight 
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The 



track was used in calculations of dynamic heights, 
results are shown in Figure 4.22 and Table XX. 

The variation of dynamic height along the track was low 
(0. 19 meters/kilometer). Across the eddy the height change 
was 0.1 meters per 150 kilometers. This height signal would 
be too small for the SEASAT altimeter to detect accurately. 
However, future generations of satellite altimeters should 
be able to detect the magnitude of the signal from this 
edd y . 

The temperature gradient in the thermocline can be an 
important factor in a model of mixed layer depth. This 
gradient cannot be predicted accurately from sea surface 
temperature observations. Therefore, a correlation between 
thermocline gradient and dynamic height was investigated. 
The gradient from one hundred meters showed the strongest 
correlation with dynamic height. The possibility exists, 
then, to predict the gradient from measured values of 
dynamic heights along the track. 

F. OVERALL OBSERVATIONS 

1 . An Edd^ on the Southeastern Part of the Track 

A warm eddy formed in the observation area. It 
began as an intrusion of warm water (the meander imaged by 
satellite on 17 November) and eventually appeared to form a 
closed circulation. Based on the November image, the eddy 
appeared to have a 90-nautical-mile (165 kilometer) major 
axis and a 60-nautical-mile (110 kilometer) minor axis. By 
the end of the experiment, the major axis had decreased to 
about 60 nautical miles, while the minor axis did not 
decrease significantly. 

The sea surface temperature and mixed layer depth 
were highly correlated across the eddy. The southeast 
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o\ 

© 



Surface Heights (dynamic meters) 

oo h- vo wn 

O O O O 




gure 4.22 



Dynamic Heights Along the Center Track. 
(17 Nov) 0/350. 
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"wall" of this eddy displayed sharp gradients of both sea 
surface temperature and mixed layer depth. The correlation 
was consistently strong enough to infer that the influence 
of the eddy was confined to the vertical area directly below 
the surface thermal manifis tation, that is, the front which 
existed between the eddy water and the surrounding water had 
a slope which approached ninety degrees, especially to the 
southeast of the eddy. 

2. Gradients in the Thermocline 



The gradient computed from the first five meters of 
the thermocline showed a negative correlation with position 
in seven of eight flight tracks. The gradient was strongest 
to the east of station nine. Station nine lies near the 
boundary between the Central Subarctic Domain and the 
Transitional Domain. 

The gradient computed from five meters was more 
correlated with position in November than in December. 
Table XXI lists the correlation coefficients for the 
gradient from five meters. 



Of the gradients computed, the gradient from five 
meters and the gradient from one hundred meters generally 
had the best correlation with position. Attempts to model 
mixed layer depth were generally more successful using the 
gradient computed from 100 meters. This is to be expected 
as the gradient computed from 100 meters tends to smooth 
errors introduced by an arbitrary starting poin 4 - (the 
strength of the five meter gradient could vary by as much as 
twenty percent depending on what depth was selected as the 
top of the thermocline. Table XXII lists the two gradients, 
by flight track, which exhibited the largest magnitude 
correlation with position on the flight track. 
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TABLE XXI 



Correlation Coefficient (r) for Gradient from 5 meters 

vs. Position 



FLIGHT 


DATE 


CORRELATION COEFFICIENT (r) 


1 




15 


No v 


-0.818 


2 




17 


No v 


-0. 583 


3 




19 


No v 


-0. 766 


4 




1 


Dec 


-0. 159 


5 




3 


Dec 


-0.478 


5 


(south) 


3 


Dec 


+0.012 


6 


5 


Dec 


-0.266 


6 


(north) 


5 


Dec 


-0.129 










TA3LE XXII 


Gradients 


Most 


Co rrel 


ated with Position (by Fligh 



FLIGHT 


DATE 


GRADIENT 


(r) 


1 


15 Nov 


5 


-0.818 






15 


-0.412 


2 


17 Nov 


5 


-0.583 






100 


+0.575 


3 


19 Nov 


5 


-0.766 






50 


+0.282 


4 


1 Dec 


5 


-0. 1 59 






50 


+ 0. 1 15 


5 


3 Dec 


100 


+0.342 






5 


-0.478 


5 (S) 


3 Dec 


100 


+0.919 






50 


+0.359 


6 


5 Dec 


1 00 


+0.781 






50 


+0.542 


6 (N) 


5 Dec 


100 


+0.524 






5 


-0. 1 29 


Not e : 









1. Track five south indicated by 5(S) and 
track six north indicated by 6 (N) . 
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3 • Se a Surface T em p er a ture 

The sea surface temperature correlated with position 
on all flight tracks. This is due in part to the system of 
specifying positions on the line. It is still interesting 
to note that the sea surface temperature increased almost 
linearly along the track until reaching the eddy. There, a 
second order fit with position was more appropriate. Table 
XXIII lists the correlation coefficient computed for sea 
surface temperature vs. position for each track. 



TABLE XXIII 

Correlation Coefficient (r) for SST vs. Track Position. 



FLIGHT DATS 



ENTIRE TRACK EDDY LOCATION 



1 




15 


Nov 


0. 3 56 


2 




17 


Nov 


0.9 09 


3 




19 


N o v 


0.569 


4 




1 


Dec 


0. 9 09 


5 




3 


Dec 


0. 8 94 


5 


(S) 


3 


Dec 


0. 9 82 


6 


5 


Dec 


0. 9 17 


6 


(N) 


5 


Dec 


0.9 54 



-0.125 

-0.416 

-0.480 

-0.648 

-0.476 

0.764 

-0.795 

0.491 



Not e: 

1. Eddy location is the along track correlation 
between buoys that would define the eddy 
boundaries (stations 4 to 24 on 17 November) . 



The objective of this study was not the formulation 
of a predictor for sea surface temperature. Observed sea 
surface temperature patterns did show when deviations from 
the expected climatology occurred. The sea surface tempera- 
ture gradients, in particular, provided very good informa- 
tion for defining frontal conditions at the eddy. Although 
values of sea surface temperature did not imply 
corresponding values of mixed layer depth with any 
reliability, the sharp gradients of surface temperature 
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similar 



along the track did occur in conjunction with 
gradients of mixed layer depth. 

Satellite-derived sea surface temperature gradients, 
then, can be used to modify sound velocity profiles when 
water mass discontinuities are observed. As a first guess, 
a single bathythermograph could provide temperatures, while 
sat el lit e- derived surface temperatures could be used to 
modify the measured vertical structure at another location. 
The sea surface temperature gradients would be useful as a 
decision aid for predicting the changes occurring in the 
mixed layer depth. A vertical temperature structure along 
the proposed track could be worked out in this manner. 



TABLE XXIV 

Correlation Coefficient (r) - MLD vs. Track Position. 



FLIGHT 


DATE 


ENTIRE 




TRACK 


EDDY 


LOCATION 


1 




15 


Nov 


- 0 . 


1 


49 


-0 


.231 


2 




17 


Nov 


- 0 . 


1 


64 


-0 


.607 


3 




19 


Nov 


- 0 . 


6 


41 


-0 


. 544 


4 




1 


Dec 


- 0 . 


6 


55 


-0 


.931 


5 




3 


Dec 


- 0 . 


5 


69 


-0 


.018 


5 


(S) 


3 


Dec 


- 0 . 


8 


45 


-0 


.527 


6 


5 


Dec 


- 0 . 


9 


27 


-0 


.954 


6 


(N) 


5 


Dec 


- 0 . 


7 


87 


-0 


. 972 


Not e 


: 


Eddy locati 


on is th 


3 


same 


as in t 


he previ 



4 . Mi xed Laver Depth 

The mixed layer depth varied along each track. 
Sometimes it would correlate with position, but seldom would 
the correlation be as high as for sea surface temperature. 
The correlation was always negative, that is the mixed layer 
shallowed toward the coast. This is to be expected from the 
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climatoioqical description of the ar<=a. The Transitional 
Domain is characterized by a weakening of the halocline r 
warmer water, and shallower isothermal surface layers. 

The sea surface temperature indicated the presence 
of a warm meander or eddy. The satellite image from 17 
November shows the thermal signature of warm water entrained 
in a cold area. Even so, the sea surface temperature 
generally did not show as strong a correlation with position 
in the eddy as did mixed layer depth. Table XXIV shows the 
correlation coefficients for mixed layer depth vs. track 
position, both for the entire track and the eddy only. This 
may provide a means of estimating mixed layer depth in warm 
eddies. Satellites can be used to define the surface 
boundaries of a warm eddy. Then mixed layer depth across 
the eddy could be predicted as a function of position in the 
edd y. 



TABLE XXV 



Correlation Coefficients (r) 



for MLD vs. S ST. 



Flight 



r (entire track) 



(eddy) 



1 




0.334 


0.970 


2 




0.170 


0.858 


3 




-0.101 


0.509 


4 




-0.418 


0.501 


5 




-0.382 


0.881 


5 


(S) 


-0.838 


-0.562 


6 


-0.546 


0.399 


6 


(N) 


-0.608 


-0.433 


Not 


e : 


The eddv column refe 
gathered only from s 
include the eddy. 


rs to corr 
tations al 



the 



of data 
crack which 
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5 



• jli xsd Layer Depth and Sea Su rf ace 'Temperature 
Cor relati ons 

Mixed layer depth and sea surface temperature showed 
little correlation in November, but the December flight 
tracks showed a strong increase in the correlation 

coefficient for the entire track. The eddy displayed a 

strong mixed- la yer-de pth se a-sur face- temperature correlation 
in November, decreasing somewhat in December. Table XXV 
summarizes the correlation coefficients for mixed layer 
depth and sea surface temperature for each day of 
observe t ions. 
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V. CONCLUSIONS AND RSCOM MSN DAT 10 NS 
A. CONCLUSIONS 

Satellites can be used to map sea surface conditions, 
particularly sea surface temperature. For ASTREX data, 
knowledge of sea surface temperature and position along the 
flight track was used no model mixed layer depth. This 
model correlated poorly with in situ observations. When 
gradient in the thermocline was added to the model, the 
correlation increased significantly. 

The gradient in the thermocline may be predicted (with 
low confidence) from satellite altimeter measurements of the 
ocean surface. A useful relationship between thermocline 
gradient and the dynamic topography variations may exist. 

The most striking result of this study was the strong 
correlation observed between sea surface temperature and 
mixed layer depth across a warm eddy. Satellite imagery 
provided the location of the eddy boundaries with a high 
degree of accuracy. The strength of this surface thermal 
front is an indication of mixed layer depth variability in 
the area. Qualitatively, a large thermal gradient at the 
surface indicates a correspondingly large change in mixed 
layer depth. Across the southeastern boundary of the warm 
eddy, the sea surface temperature and mixed layer depth 
changed sigini f icant ly in the same sense (both increased 
together when tracked from the northwest to the southeast). 

Uncorrected satellite temperatures in the area displayed 
the same variations along the track as did observed AXBT sea 
surface ■ t em pe nature s. The variations of sea surface temper- 
ature (thermal gradients at the surface) contain the infor- 
mation needed to locate anomalies such as the warm eddy 
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observed in the southeastern part of the flight tracks. The 
NOAA-6/7 satellites can resolve 0.5 degree Kelvin tempera- 
tures at one-kilometer resolution. Frontal areas can be 
identified and accurately located. The AXBT's provided 
temperature values at 55-kilometer intervals, which were 
sufficient to observe that the eddy did exist, but were not 
sufficient to pinpoint the eddy boundaries. Satellite 
imagery, providing temperature values at one-kilometer 
intervals, can be used no specify the boundaries of a front 
or eddy more accurately. 



B. RECOMMENDATIONS 



In future experiments in the ASTREX region, 
consideration should be given to investigating satellite 
located- anomalies such as the warm eddy observed during 
ASTREX. Until a satellite altimeter becomes available, 
dynamic heights must be calculated indirectly from other 
observed variables. Therefore, the salinity should be 
measured in any future test area. 



Various 


models of mi 


xed layer depth 


must be explored. 


In addition 


to modeling 


the mixed layer 


as a function of 


position , 


sea surface 


temperature and 


gradient in the 


then mod in e 


, additional 


factors should 


be added. Thes* 



factors should include wind mixing in the mixed layer based 
on observed wind speeds in the area, dynamic heights 
associated with tracks examined, and possibly observed sea 



state in the area. 

Another model which could prove v<=ry useful would 
include the effects of wind speed and duration on the mixed 
layer depth. Wind was not used as a factor in this study, 
but it is an important element in the formation cf the mixed 
layer. Between the November and December flights, wind 
speed increased and so did mixed layer depth. Instruments 
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such as 
for an 
function 
tern perat 
of mixed 



the SMMR on NIMBU3-7 can provide scaler 
area. The formulation of mixed layer 
of magnitude and duration of wind, s 
ures, and other factors may provide the 
layer depth. 



wind 



depth 



ea s 
best. 



fields 
as a 
urf ace 
model 
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APPENDIX A 

VARIAT ION OF THERKOCLr N F GR ADI ENTS WITH POSITION FOR 

NOVEMBER 

— 

Gradient 

(degrees Celsius/meter) 

0. 560 




0. 16 C 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



» 1 



Figure A- 1 Gradient fro* 5 meters vs. Position (15 Nov). 
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r 



Gradient 

(deg rees Celsius/meter) 
0. 26 0 




0 . 160 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Milas from Origin) 



Figure A. 2 



Gradient fro® 15 meters vs. Position (15 Nov). 
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r 



T 



Gradient 

(degrees Celsius/m eter) 
0 . 1120 




0. 0720 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure A. 3 Gradient from 50 meters vs. Position (15 Nov). 
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Gradient 

(degrees Celsius/meter) 




0. * ” 1 50 . 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure A. 4 Gradient froa 100 meters vs. Position (15 Nov). 
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r 



Gradient 

(degrees Celsius/meter) 




Figure A. 5 Gradient fro* 5 meters vs. Position (17 Nov). 
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Gradient 

(degrees Celsius/meter) 
0. 270 




0 . 120 

+ + + +- 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure A. 6 Gradient 



15 meters vs. Position (17 Nov). 
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Gradient 

(degrees Celsius/meter) 




+ _+ + + +- 

0 . 150 . 300 . 450 . 600 . 

TRACK POSITION (Nautical Miles from Origin) 



Figure A. 7 



Gradient from 50 meters vs. Position (17 Nov). 



Gradient 

(degrees Celsius/meter) 
0.0700 




0. 0350 

+ + + + +- 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Pigure A. 3 Gradient from 



I v/u 
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Gradient 

(decrrees Celsius/m eter) 
0. 560 




Figure A. 9 



Gradient froi 5 meters vs. Position (19 Nov). 
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Gradient 

(degrees Celsius/meter) 
0. 260 




0. 160 

0. 100. 200. 300. 400. 

TRACK POSITION (Nautical Milas from Origin) 




Figure A. 10 



Gradient from 15 meters vs. Position (19 Nov). 



Gradient 

(degrees Celsius/meter) 

0 . 1120 




0. 0720 

+ + ♦ + + - 

0. 100. 200. 300. 400. 

TRACK POSITION (Nautical Miles from Origin) 



Pig are A. 11 



Gradient from 50 meters vs. Position 



(19 Nov). 



r 



Gradient 

(degrees Celsius/meter) 




0. 100. 200. 300. 400. 

TRACK POSITION (Nautical Miles from Origin) 



Figure A. 12 Gradient from 100 meters vs. Position (19 Nov). 
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APPENDIX B 



VARIATION OF THERMOCLINE GRADIENTS WITH POSITION FOR 

DECEMBER 

i 1 



Gradient 

(degrees Celsius /meter) 
0.560 




0. 160 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



— 



Figure B. 13 Gradient fro« 5 meters vs, Position (1 Dec). 
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Figure B. 1 4 



Gradient from 15 meters vs. Position 



(1 Dec) 



Gradient 

(degrees Celsius/meter) 

0 . 1120 




0. 0720 

+ + + + + - 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



i 1 



Figure B. 15 



Gradient from 50 meters ts. Position 



(1 Dec) 



1 



Gradient 




o7~ 150 7" 300. 4 5 o7 600. 

TRACK POSITION (Nautical Miles from Origin) 



« 



Figure B.16 



Gradient froa 100 meters vs. Position 



(1 Dec). 



r 



Gradient 

(degrees Celsius/meter) 
0.50 




Figure B. 17 



Gradient fro* 5 meters vs. Position 



(3 Dec) 



Gradient 




TRACK POSITION (Nautical Miles from Origin) 



Figure B. 18 Gradient froa 15 meters vs. Position (3 Dec) 



r 



Gradient 

(degrees Celsius/m eter) 



0. 100 




+ + + + + - 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure B. 19 



Gradient fro* 50 meters vs. Position 



(3 Dec). 



r 



Gradient 

(degrees Celsius/meter) 

0.060+ 

0. 050 

0.040~ 

0. 030 

0 . 020 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure B.20 Gradient from 100 meters vs. Position (3 Dec). 
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Gradient 

(degrees Celsius/meter) 




Figure B.21 



Gradient fro* 5 meters vs Position (3 Dec 
South). 
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Gradient 

(degrees Celsius/meter) 




0. 090 

0. 150. 300. U50. 600. 

TRACK POSITION (Naurical Miles from Origin) 



Figure B.22 Gradient from 15 meters vs Position (3 Dec 
South). 
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Gradient 

(degrees Celsius/meter) 
0.09 0 




0.040 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure B.23 



Gradient from 50 meters vs Position (3 Dec 
South)* 
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Gradient 

(degrees Celsius/meter) 
0.0560 




0. 0210 

+ + + + + - 

0. 150. 300. 450. 600. 

TBACK POSITION (Nautical Miles from Origin) 



Figure B.24 Gradient froi 100 meters vs Position (3 Dec 
- South). 
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Gradient 

(degrees Celsius/m etsr) 
0. 400 




0. 150 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 

I i 

Figure B. 25 Gradient fro* 5 meters vs. Position (5 Dec). 
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Gradient 

(degrees Celsius/meter) 
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+ + + + +- 
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TRACK POSITION (Nautical Miles from Origin) 



Figure B. 26 



Gradient from 15 meters vs. Position (5 Dec), 
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Figure B. 27 Gradient from 50 meters vs. Position (5 Dec), 
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r 



Gradient 

{degrees Celsius/meter) 



0. 060 




Pigure B.28 



Gradient from 100 meters vs. Position (5 Dec). 
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Gradient 

(degrees Celsius/meter) 
0. 490 
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+• — — — — t- — — — + —————— — — + — — — 

0. 150. 300. 450. 600. 

TRACK POSITION (Nautical Miles from Origin) 



Figure B.29 



Gradient froa 5 meters vs. Position (5 Dec 
North). 
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Gradient 

(degrees Celsius/meter) 
0. 240 




+ + + + + 



0. 150. 300. 450. 600. 
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Figure B.30 



Gradient from 15 meters vs. Position (5 Dec - 
North). 



130 



Gradient 
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Figure B.31 



Gradient from 50 meters vs. Position (5 Dec - 
North). 
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Figure B.32 



Gradient from 100 meters vs. Position (5 Dec 
North). ' 
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APPENDIX C 



TABLES OF NORMALIZED GRADIENTS IN THE THERMOCLINE 



TABLE XXVI 





Normalized Grad 


ients for 


Flight Track 


5 meter 


15 meter 


50 meter 


100 meter 


0.3980 


0.1920 


0. 0818 


0.0479 


0.5340 


0.2440 


0. 0948 


0.0549 


0.4 96 0 


0.2027 


0. 0802 


0.0530 


0.3880 


0.2220 


0. 0898 


0.0566 


0 .5040 


0.2240 


0. 0844 


0.0552 


0.3400 


0.2427 


0. 1054 


0.0649 


0.3120 


0.1720 


0. 0802 


0.0506 


0.2260 


0.2093 


0. 1070 


0.0622 


0 .2820 


0.1793 


0. 0924 


0.0592 


0.2300 


0.1907 


0. 0780 


0.0465 


)te: 


Gradients were com put 


ed from z 


he temperature 


change 

column 


over the distance speci 
head! ng. 


fied in the 
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TABLE XXVII 



Normalized Gradients for Plight Track Iwo. 



5 meter 


15 meter 


50 meter 


100 met 


0.4940 


0.2060 


0. 0780 


0.0425 


0.3700 


0. 1753 


0. 0774 


0.0423 


0.3420 


0.1787 


0. 07 24 


0.0424 


0 .5 160 


0.2260 


0. 0946 


0.0519 


0 .5 160 


0.2240 


0. 0894 


0.0501 


0.3340 


0.2467 


0. 0952 


0.0552 


0 .4 160 


0.2480 


0. 0964 


0.0525 


0. 2380 


0.1973 


0. 0872 


0.0495 


0.3500 


0. 1920 


0. 0780 


0.0495 


0.5300 


0.2233 


0. 0928 


0 . 0530 


0 .2860 


0.1920 


0. 0814 


0.0564 


0 .3380 


0.2167 


0. 0866 


0.0570 


0.5320 


0.2420 


0. 1012 


0.0622 


0.3280 


0.1953 


0. 0818 


0.0533 


0.2320 


0.2367 


0. 1040 


0.0624 


0.3020 


0.2313 


0. 09 90 


0.0611 


0.3660 


0.2560 


0. 0968 


0.0601 


0. 2960 


0.2473 


0. 10 02 


0.0604 


0.3600 


0.2333 


0. 1140 


0.0670 


0 .3020 


0.1880 


0. 0954 


0.0653 


0. 1 860 


0.1760 


0. 0932 


0.0630 


0 . 2660 


0. 1567 


0. 0680 


0 . 0397 


0.2160 


0.1300 


0. 0856 


0.0535 



TABLE XXVIII 

Normalized Gradients for Flight Track Three. 



5 meter 


15 meter 


50 meter 


100 mete 


0 .4960 


0.2027 


0. 0802 


0.0530 


0 .5 34 0 


0 . 2440 


0. 0948 


0.0549 


0.3880 


0.2220 


0. 0898 


0.0566 


0.3980 


0. 1920 


0. 0818 


0.0479 


0.5040 


0.2240 


0 . 0844 


0.0552 


0 .3400 


0.2427 


0. 1054 


0 . 0649 


0.2300 


0.1907 


0. 0780 


0.0465 


0.3120 


0.1720 


0. 0802 


0.0506 


0.2820 


0. 1793 


0. 09 24 


0.0592 


0 .2260 


0.2093 


0. 1070 


0 . 0622 


0.2260 


0.2093 


0. 1070 


0.0622 


0 .2820 


0.1793 


0. 0924 


0.0592 


0.3400 


0.2427 


0. 1054 


0.0649 


0.2300 


0.1907 


0. 0780 


0.0465 
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TABLE XXIX 



Normalized Gradients for Flight Track Fou 



5 meter 


15 meter 


50 meter 


100 meter 


0 .5340 


0.2440 


0 .0948 


0.0549 


0.4 96 0 


0.2027 


0 .0802 


0. 0530 


0.3980 


0. 1920 


0 .081 8 


0.0479 


0.5040 


0.2240 


0 . 0844 


0.0552 


0 . 2300 


0.1907 


0 .0780 


0. 0465 


0.3880 


0.2220 


0 . 0898 


0. 0566 


0 .5040 


0.2240 


0 .0844 


0.0552 


0.2820 


0. 1793 


0.0924 


0. 0592 


0.3880 


0.2220 


0.0898 


0.0566 


0.2260 


0.2093 


0 . 1070 


0.0622 


0.3400 


0.2427 


0 .1054 


0.0649 


0 .3 120 


0. 1720 


0.0802 


0. 0506 


0 . 3 120 


0. 1720 


0.0802 


0. 0506 


0.3400 


0.2427 


0. 1054 


0.0649 


0. 2260 


0.2093 


0 . 1070 


0.0622 


0 .3980 


0. 1920 


0.081 8 


0. 0479 


0. 2820 


0. 1793 


0 .0924 


0.0592 


0.5340 


0.2440 


0.0948 


0.0549 


0 .496 0 


0.2027 


0.0802 


0.0530 


0.5340 


0.2440 


0.0948 


0. 0549 


0.2300 


0. 1907 


0 .0780 


0. 0465 



TABLE XXX 

Normalized Gradients for Flight Track Five 



5 meter 


15 meter 


50 meter 


100 meter 


0 .2840 


0.1253 


0. 0440 


0.0249 


0 . 1 560 


0.0953 


0. 0448 


0. 0243 


0.4560 


0. 1807 


0. 0660 


0.0382 


0.3600 


0. 1507 


0. 0648 


0.0337 


0 .3060 


0. 1773 


0. 0742 


0.0436 


0.1620 


0. 1380 


0. 0602 


0.0414 


0.2840 


0.1627 


0 . 0840 


0.0520 


0 .2820 


0. 1447 


0. 0670 


0. 0451 


0.2700 


0.1820 


0. 0974 


0.0546 


0.2100 


0.1593 


0. 0934 


0.0597 


0.0080 


0. 0853 


0. 0716 


0. 0448 



(Center) 
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rABLE XXXI 



Normalized Gradients for Flight Track Five 



(South) . 



5 me+er 

0 .2460 
0.2860 
0.3380 
0.2000 
0.1840 
0 .3580 
0.3240 
0.2140 



15 meter 

0. 1707 
0.1133 
0.1253 
0. 1687 
0.1733 
0.2407 
0. 1707 
0.1400 



50 meter 

0. 0646 
0. 0480 
0. 0616 
0. 0668 
0. 0696 
0. 0840 
0. 0840 
0. 0818 



100 meter 

0.0342 
0.0262 
0.0355 
0. 046 1 
0. 0489 
0.0476 
0.0512 
0.0522 



TABLE XXXII 

Normalized Gradients for Flight Track Six 



(Center) . 



5 meter 


15 meter 


0.2920 


0.1220 


0. 3200 


0.1507 


0.3520 


0.1713 


0.2220 


0.1640 


0 . 2780 


0.1633 


0.2200 


0.1453 


0.3540 


0.2033 


0.3580 


0.2033 


0. 1780 


0. 1400 


0.3540 


0.2113 


0. 1680 


0.0780 



50 meter 


100 meter 


0. 0404 


0.0232 


0. 0582 


0.0283 


0. 0694 


0.0375 


0. 0684 


0.0342 


0 . 0734 


0.0425 


0. 0680 


0.0430 


0. 0800 


0.0483 


0. 0920 


0.0544 


0. 0850 


0.0588 


0. 0758 


0. 0440 


0. 0638 


0.0421 
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TABLE XXXIII 



Normalized Gradients for Flight Track Six 



5 meter 


15 meter 


0 . 2440 


0.1060 


0 .2720 


0.1547 


0 . 4 06 0 


0.1787 


0.4480 


0. 1927 


0.1820 


0. 1440 


0.3400 


0.1540 


0.1900 


0.1313 


0.1840 


0. 1520 


0.3280 


0.2280 


0 .3720 


0.2007 


0 . 2 24 0 


0 .0960 



50 meter 


100 meter 


0. 0428 


0.0237 


0. 0678 


0.0369 


0. 0730 


0.0377 


0. 0758 


0.0423 


0. 0478 


0.0269 


0. 06 56 


0.0449 


0. 0606 


0.0352 


0. 0686 


0.0489 


0. 0792 


0.0462 


0. 0774 


0.0471 


0. 0456 


0.0359 



(North) . 
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APPENDIX D 



MODEL EQUATIONS FOR EACH FLIGHT 



TABLE XXXIV 



Model Equations for MLD (15 November) 



Model Variables Equations 



SST 




MLD 


SST , 


Posit ion 


MLD 


SST. 


Posit ion. 


MLD = 



Gradient (100) 



- 1 22 

- 0.0918 (Position) 
+ 14.6 (SST) 

- 1 05 

- 0 .0718 (Position) 
♦ 10.5 (SST) 

+ 62.3 (GRAD ( 100) ) 



Not es z 

1. SST is sea surface temperature. 

2. Position refers to station position 
relative to the most seaward station. 

3. GRAD (100) refers to the gradient 
computed from the upper 100 meters of 
the thermocline. 
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TABLE XXXV 



Model Equations for MLD (17 November) 



Model Variables 
SST 

SST , Position 



SST, Position 
Gradient (15) 



9 



SST, Position. 
Gradient (100) 



Equations 



MLD 

MLD 



+ 4 3. 7 + 1.01 (SST) 



57.7 

0.0865 {Position) 
10.9 (SST) 



MLD = + 6 2. 8 

- 0 .098 (Position) 
♦ 1 2. 3 (SST) 

- 4 3. 1 (GRAD (15) ) 



MLD = - 8 8.2 

- 0 . 1 15 (Position) 
+ 17.3 (SST) 

+ 8 18 (GRAD (100) ) 



TA3LE XXXVI 

Model Equations for MLD (19 November) 



Model Variables 


Equations 




SST 


MLD = + 


6 8.7 


- 0.86 (SST) 


SST, Position 


MLD = + 
+ 


17.9 

0 .055 1 (Position) 
3.71 (SST) 


SST, Position, 
Gradient (5) 


MLD = - 

+ 

+ 


7.46 
0 .045 
4 .58 
3 1.2 


(Position) 

(SST) 

(GRAD (5) ) 


SST, Position. 
Gradient (100) 


MLD = + 


1 8.2 
0 .055 
3 .72 
8.62 


(Position) 

(SST) 

(GRAD ( 100) ) 
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TABLE XXXVII 



Model Equations for MLD (1 December) . 



Model Variables 
SST 

SST , Position 

SST , Position, 
Gradient (50) 



SST, Position, 
Gradient (100) 



Equations 



MLD = + 92. 5 -2.03 (SST) 

MLD = + 2 8.9 

- 0.061 (Position) 
+ 4.93 (SST) 



MLD = + 11.3 

- 0 .06 (Position) 
+ 4.62 (SST) 

+ 2 35. (GRAD (50) ) 



MLD = + 16.3 

- 0.057 (Position) 
+ 4.28 (SST) 

+ 3 44 (GRAD (100) ) 



TABLE XXXVIII 

Model Equations for MLD (3 December, Center) 



Model Variables 


Equations 


SST 


MLD = + 


9 5. 5 - 2. 48 (SST) 


SST, Position 


MLD = + 
+ 


4 0.61 

0 .0683 (Posi- 1 - ion) 
4.10 (SST) 


SST, Position, 
Gradient (15) 


MLD = + 

+ 

+ 


5 2. 0 

0 .043 (Posi ti on) 
0.95 (SST) 

118.2 (GRAD (15)) 


SST, Position, 
Gradient (100) 


MLD = + 
+ 


7 1.7 

0 . 064 (Posi tion) 
3.61 (SST) 

1 384 (GRAD (10 0) ) 
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TABLE XXXIX 



Model Equations for MLD (3 December, South). 



Model Variables 


Equations 




SST 


MLD = + 


95.5 


- 2.48 (SST) 


SST, Position 


MLD = + 
+ 


4 0.61 

0 .0683 (Position) 
4.10 (SST) 


SST , . Posit ^on. 
Gradient (50) 


MLD = + 
+ 


1 27 
0 .021 
0 .85 
9 36.6 


(Position) 

(SST) 

(GRAD (50) ) 


SST. Position. 
Gradient (100) 


MLD = + 
+ 


1 02 
0 .044 
1 .530 
1 03.0 


(Position) 

(SST) 

(GRAD (100) ) 



TAELE XL 

Model Equations for MLD (5 December, Center). 



Model Variables 
SST 

SST , Position 

SST , Position, 
Gradient (5) 



SST, Position. 
Gradient (100) 



Equations 

MLD =+111- 

MLD = + 12.40 
- 0.103 
+ 8.18 



MLD = + 3 .23 
- 0.10 
+ 7.92 
+ 20.9 



MLD = + 20.2 
- 0.099 
+ 6.70 
+ 187.4 



3.38 (SST) 



(Position) 

(SST) 



i Posit ion) 
’GRAD ( 5) ) 



(Position) 

(SST) 

(GRAD (100)) 
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TABLE XLI 



Model Equations for MLD (5 December, North) 



Model Variables 


Equations 




SST 


MLD = + 


1 23 - 


5.01 (SST) 


SST , Position 


MLD = - 
+ 


4 3.7 
0 .12 
1 3. 1 


(Posit io n) 
(SST) 


SST, Position, 
Gradient (5) 


MLD = + 

+ 

+ 


9.23 

0.10 i 

7 .92 
20.9 


(Position) 

(SST) 

(GRAD (5) ) 


SST, Position. 
Gradient ( 1 00 f 


MLD = - 

+ 

+ 


3 4.0 

0 . 1 14 

1 1.42 
1 74.7 


(Position) 

(SST) 

(GRAD (100) ) 
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